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Theory  of  time-dependent  propagation  in  multimode  lightguides* 

Bruno  Crosignani*  and  Charles  H.  Papas 
California  Institute  of  Technology,  Pasadtna,  California  91I2S 
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A  coupled  system  of  equations  governing  the  propagation  of  a  signal  in  a  statistical  ensemble  of  multimode 
optical  fibers  is  presented.  It  describes,  besides  the  usual  average  modal  powers,  the  evolution  of  the 
interference  terms  between  the  mode  amplitudes  and  of  the  modal  power  fluctuations.  Our  procedure  allows 
us  to  treat  the  general  nonstationary  nonmonochromatic  case  of  an  arbitrary  signal  fed  into  the  lightguide  by 
a  source  possessing  a  finite  spectral  bandwidth.  The  introduction  of  modal  power  I’.ucaw'ions  permits  ua  to 
establish  a  theorem  connecting  the  value  of  the  modal  power,  averaged  over  the  fiber  ensemble,  with  the 
actual  one  concerning  a  single  fiber.  These  two  values  coincide,  in  the  polvc*--  ~  v-  case.  <br  large  values 
of  the  fiber  length,  thus  providing  the  main  result  of  the  paper,  that  is  the  justification  of  the  statistical 
approach  to  the  problem  of  propagation.  Furthermore,  the  analysts  of  the  interference  terms  presents 
evidence  for  the  difference  between  the  propagation  of  an  amplitude-modulated  and  a  frequency-modulated 
signal. 


I.  INTRODUCTION 

Multimode  optical  fibers  with  large-diameter  core  are 
particularly  useful  in  connection  with  sources  possessing 
intermediate  degrees  of  spatial  coherence,  such  as  light- 
emitting  diodes  and  laser  diodes.  These  sources  have  a 
finite  bandwidth  6v,  and,  accordingly,  the  problem  of 
propagation  in  the  fiber  cannot  be  simplified  by  assum¬ 
ing  a  monochromatic  input.  Indeed,  a  theory  which  takes 
into  account  the  lack  of  monochromaticity  (or  the  finite 
coherence  time  t,  =  1/6  v)  is  necessary. 

The  description  of  the  propagation  of  an  electromag¬ 
netic  signal  in  a  multimode  lightguide  in  the  presence  of 
mode  coupling  is  a  statistical  one.  This  is  a  consequence 
of  the  fact  that  the  guiding  structure  unavoidably  pre¬ 
sents  some  imperfections  responsible  for  the  coupling 
among  the  modes,  whose  randomness,  together  with  the 
analytical  complexity  of  the  exact  deterministic  problem, 
leads  in  a  natural  way  to  the  statistical  approach.  It 
consists  of  an  introduction  of  an  ensemble  of  macroscop- 
ically  similar  fibers,  which  differ  among  themselves  only 
for  the  random  behavior  of  the  microscopic  imperfec¬ 
tions,  and  we  consider  the  values  of  ail  the  significant 
quantities  averaged  over  this  ensemble. 1  In  this  frame, 
the  main  emphasis  has  been  placed  on  the  evaluation  of 
the  average  power  carried  by  each  propagating  mode, 
both  in  the  case  of  a  monochromatic1  and  of  a  poly¬ 
chromatic  signal.  *’* 

However,  it  is  hardly  necessary  to  remark,  as  al¬ 
ways  when  dealing  with  a  statistical  approach,  about  the 
relevance  of  the  fluctuations  over  the  average  value, 
since  only  their  (relative)  smallness  assures  that  the 
average  value  of  a  given  quantity  can  be  safely  assumed 
as  the  actual  one.  The  fluctuations  have  been  evaluated, 
only  in  the  monochromatic  case, 1,4  with  a  rather  disap¬ 
pointing  result.  More  precisly,  it  has  been  found  that 
the  normalized  mean-square  value  ot  the  power  carried 
by  each  mode  approaches  2  for  long-traveled  distances, 
so  that  there  is  a  100%  uncertainty  on  the  actual  value  of 
this  quantity  in  a  single  fiber  once  the  average  value  is 
known. 

The  above  discussion  gives  evidence  to  the  necessity 


of  refining  the  theory  of  propagation  in  a  multimode  opti-  * 
cal  fiber  in  two  ways.  First,  the  existing  treatment 
only  deals  with  the  average  power  per  mode,  which  is 
connected  with  the  square  value  of  the  corresponding 
mode  amplitude,  while  a  more  complete  description  al¬ 
so  requires  the  study  of  the  evolution  of  the  correlation 
terms  between  different  mode  amplitudes.  In  the  course 
of  the  paper  it  wilt  be  shown  that  these  are  the  only 
terms  influenced  by  a  frequency  modulation  of  the  input 
signal,  an  operation  which,  conversely,  does  not  affect 
the  behavior  of  the  power  per  mode.  Second,  one  needs 
a  more  complete  investigation  of  the  fluctuations,  cap¬ 
able  of  covering  the  nonmonochromatic  case.  In  fact, 
it  will  be  proved  that  the  finite  bandwidth  of  the  source 
has  the  positive  consequence  of  reducing  the  magnitude 
o!  the  fluctuations  so  that,  for  long-traveled  distances, 
the  average  power  per  mode  does  coincide  with  the  actual 
value  in  a  single  fiber. 

The  mentioned  results  are  obtained  in  the  frame  of  a 
general  theory  of  nonstationary  propagation,  dealing  with 
a  carrier  possessing  a  finite  bandwidth  and  modulated 
in  an  arbitrary  way,  in  which  the  complete  description 
of  the  propagation  and  of  the  associated  fluctuations  Is 
obtained  in  the  form  of  a  closed  system  of  differential 
equations  connecting  all  the  correlation  terms  between 
mode  amplitudes  up  to  the  fourth  order.  The  present 
treatment  generalizes  that  of  Marcuse, 1  since  it  is  ap¬ 
plicable  to  the  propagation  of  a  polychromatic  signal,  a 
case  which  substantially  differs  from  that  relative  to  the 
monochromatic  signal,  as  it  is  for  example  shown  by  the 
drastic  change  in  the  behavior  of  the  fluctuations,  which 
Improves  the  predictability  of  the  actual  value  ol  the 
single-mode  power  in  the  case  of  finite  bandwidth. 

II.  PROPAGATION  IN  A  SINGLE  FIBER 

The  transverse  part  of  a  monochromatic  electromag¬ 
netic  field  ot  angular  frequency  at,  propagating  in  the  a 
direction  along  the  axis  of  a  (nearly)  cylindrical  dielec¬ 
tric  fiber,  can  be  expressed1 

■(r,  a,  t)-£  B„(r,  w)c„(a,  w) 

M 

xexp(fot- «fl,(<»>)a] ,  (1) 
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where  B„(r,  w)  is  the  transverse  electric  field  vector  of 
the  mth  mode  of  the  ideal  cylindrical  guiding  structure, 
$»(<*>)  is  its  relative  propagation  constant,  and  the  cm’ s 
are  suitable  expansion  coefficients.  In  writing  Eq.  (1), 
we  have  assumed  that  the  electric  Held  can  be  expressed 
in  terms  of  the  forward-traveling  guided  modes  alone; 
and  we  have  allowed  the  expansion  coefficients  cm  to  be 
s  dependent  in  order  to  take  into  account  the  departure 
of  the  fiber  shape  from  the  ideal  one  to  which  the  mode 
configurations  B*  pertain.  Here  r  =  (x,  y)  indicates  the 
coordinates  in  a  plane  transverse  to  2  and  one  needs  to 
consider  only  the  transverse  part  of  the  electric  field 
E(r,  2,  t)  and  of  the  magnetic  field 

H(r,  2,  t)  =  ^  **».(*,  <*>)  c„(2,  to) 

m 

xexp[tu>t  -  (2) 

for  the  purpose  of  evaluating  the  power  carried  by  the 
field. 

Equations  (1)  and  (2)  can  be  easily  generalized  to  the 
nonmonochromatic  case.  In  particular,  if,  as  it  is 
usually  the  case,  the  bandwidth  do)  of  the  field  fulfills 
the  relation 


xexp{i(u>-w,)f-i[A.(w)-dm(u;,)]2}»23'P-(i.  0  . 

m 

(8) 

where  Pm(z,  t )  is  a  positive  quantity,  as  it  is  immediate¬ 
ly  verified  by  inspection,  which  can  be  interpreted  as  the 
electromagnetic  power  carried  by  the  mth  mode. 

If  mode  coupling  is  absent,  one  has 
cm(.z,  <»)=c„(0,  j  dle~iut 

X  JL  dx  f  rfye«‘E(r’  °»  <>xHS(r,  “0) 

*17?  C  dte',U,eIU°,jrlt) 

*  JL 1- dy  e*  ‘o(r)x  Hi(r’ w<,)  •  (9> 

having  assumed,  for  the  sake  of  simplicity,  the  electric 
field  at  2  =  0  to  be  of  the  form 

E(r,  0,  <)=e‘«o‘.F«)G(r).  (10) 

Taking  advantage  of  Eq.  (3),  one  can  approximately 
write 


6u)/u)„  « 1 , 


(3)  &.(<*>)  =  &.(«•)  +  (<*>-  «p)/  Vm  , 


(11) 


where  u>q  is  the  central  frequency,  one  can  write 
B(r,  2,  0  =  £  Em(r,  u>„)  j  ducm(z,  <*>) 

xexp(iw*-idm(a))2]  (4) 

and  an  analogous  expression  for  H.  The  complex  Poynt- 
ing  vector  is  given  by5 

S=i<BxH*>„,  (5) 

where  (  )„  indicates  the  averaging  operation  over  the 
fluctuations  of  the  source.  It  is  worth  noting  that  the 
electromagnetic  field  possesses,  in  general,  fast  fluctua¬ 
tions  associated  with  the  carrier  (that  is  with  the  source) 
and  a  slow  variation  associated  with  the  superimposed 
signal.  The  averaging  operation  indicated  as  <  >„  re¬ 
fers  to  the  statistical  ensemble  pertinent  to  the  source. 
Whenever  the  source  is  stationary,  and  the  ergodic  hy¬ 
pothesis  is  justified,  *  the  quantity  P,  in  Eq.  (6)  repre- 
sentsthe  power  averaged  over  a  time  of  the  order  of  2ir/ 
fiu>  and  still  retains  the  slow-time  variation. 

The  power  carried  by  the  field  through  an  arbitrary 
fiber  section  z  *  const  can  be  expressed 

'P'*ReJ-  dX  J-  (8) 

where  e(  is  the  unit  vector  in  the  2  ulrection  and  "Re” 
means  "real  part  of.  ”  Taking  advantage  of  the  orthogo¬ 
nality  relations  between  the  guided  modes1 

dx  «o)xH*(e,  Wpj-ipa*,,  (7) 

where  P  is  a  positive  normalization  constant  and  is 
Kronecker’s  symbol,  one  obtains  turn  ugh  Eqs.  (4)— (7), 

P<mP  T'.  j  dw  J  du‘  (cm(z,  u>)ctU, 


with 


v;MdA.(u>)/du>L.U0,  (12) 

and  Eq.  (8)  yields 


*  |  dx  dy G(r)xH*(r,  u>o).e,|  . 


(13) 


Equation  (13)  expresses  the  fact  that  the  modulation 
imposed  at  z  -  0  travels  in  each  mode  with  the  mode 
characteristic  group  velocity  V„ ;  in  particular,  that 
frequency  modulation,  which  affects  the  phase  of  F(t), 
is  ineffective  in  influencing  the  total  power  carried 
through  a  fiber  section. 


In  order  to  satisfactorily  treat  the  situation  in  which 
mode  coupling  is  present,  one  has  to  resort  to  an  en¬ 
semble-average  approach.  More  precisely,  one  intro¬ 
duces  a  statistical  ensemble  of  macroscopically  equiva¬ 
lent  fibers  possessing  microscopic  random  imperfections 
and  then  considers  the  ensemble  average,  Indicated  by 
the  symbol  (•••).  Accordingly,  all  the  significant  quanti¬ 
ties  will  be  expressed  by  means  of  two  statistically  in¬ 
dependent  averaging  operations,  the  first  of  which,  in¬ 
dicated  by  (  )m,  refers  to  the  source,  while  the  second. 
Indicated  by  angular  brackets,  is  associated  with  the  in¬ 
troduction  of  the  statistical  fiber  ensemble.  Thus  one 
has,  for  example,  the  total  average  power  carried 
through  a  fiber  section  given  by 


<P,>-]T<PJz,  t )> 

-f>£  j  dw  J  dw'«cm(2,  w)c*(2,  </)>„> 

X  exp{f(u>  -  ui'X  -  i[fl.(u>)  -  A.(ui')]  2}  .  (14) 
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Ill.  STATISTICAL  THEORY  OF  TIME-DEPENDENT  . 
PROPAGATION  IN  AN  ENSEMBLE  OF  FIBERS 

In  order  to  study  the  space-time  evolution  of  CP„(z,  1)) 
in  the  liber,  one  has  to,  according  to  Eq.  (14),  investi¬ 
gate  the  behavior  of  the  quantity  ((c„(z,  u >)c*(z,  oi')>„>. 
The  starting  point  is  the  set  of  equations  describing  the 
deterministic  evolution  of  the  c„’s  in  a  single  fiber. 
They  read1 

L  4*U,  '■>),  rn  =  1,  2 . 

(15) 

where 

A.iU,  <*>,  oi0)**-|(*,  wo) 

xexp{t[6*(ui)-6,(ui)]r} ,  (16) 

the  A^’s  being  suitable  coupling  coefficients  obeying  the 
relation 

Kml  -  -  A*.  .  (17) 

An  equation  completely  analogous  to  Eq.  (15)  holds 
obviously  true  for  oi'  and  reads 

"=  22  A./U.  to',  ui0)  c,(z,  co') ,  m=  1,  2, ...  . 

(16) 

From  Eqs.  (15)  and  (18)  one  can  immediately  derive 


the  following  set  of  equations: 

■£  [cite.  oi')e,(z,  to)] 

=  22  to,  !O0)c*U,  to')c,U,  to)  (19) 

ft 

+  ^(*.  to',  to#)  c,(z,  u>)  ef(z,  «')}  ,  m,j*  1,  2, ...  . 

At  this  point,  one  can  follow  a  procedure  analogous  to 
that  Introduced  in  the  monochromatic  case  (ui=u/)  (Ref. 
4);  that  is,  one  formally  integrates  Eq.  (19)  and  then 
substitutes  in  the  resulting  expression  for  c£(z,  w')c,(z, 
to)  the  formal  solution  for  the  quantities  c£(z,  u>')c,(z, 
io)  and  Cj(z,  u)c{(z,  o>').  After  doing  this,  one  can  per¬ 
form  the  double-averaging  operation  on  the  resulting 
equations. 

If  the  coupling  process  is  stationary  to  order  two  in  the 
z  direction,  that  is, 

<o0)  A„(z',  o>*)>  =  Fmjr.(  |  z  -  z'  | ) ,  (20) 

and  if 

(A„j(z,  !O0)>  =0,  (21) 

by  assuming  that  the  mode-field  amplitudes  do  not  ap¬ 
preciably  vary  over  distances  comparable  with  the  typi¬ 
cal  correlation  length  D  of  the  correlation  function  de¬ 
fined  in  Eq.  (20),  one  obtains 


~  A.,(z,  oi,  ui')=  £22  **rU.  ui')exp{t[8l,(ui)-&{o')3*}  <JCml(z')Klt{0)) 

x exp{-  i[6,  (oi)  -  &(ui)]  z'jrfz'  + Art(z,  oi,  oi') exp  (t[A.(u>)  -  6,(oi)  +  A,(ui')  -  6,(u/)]  z} 

x  f“(A.i(z')K:,(0)>exp{-i[A(w')-6r(a'')]z'}dz'+A,>(z,  oi,  u,') 

JO 

xeq){i[^(w')- 8,(oi')  +  A.(uj)  -  6,(oi))z}  f  (Aj/z')A^( 0)>  exp{-  t[6»(ui)  -  6»(o>)] z^dz' 

J  0 


+X—(z,  u>,  oj')  exp{i  [&(ui')  -  8,(ui')]  z 


where  we  have  defined 

X„(z,  01,  0>')»«c.(z,  ui)e?(z,  oi')>„>;  (23) 

and  the  dependence  on  ui0  of  the  coupling  coefficients 
has  been  omitted. 

It  can  now  be  observed  that  only  nonoscillatory  expo¬ 
nential  terms  in  z  contribute  appreciably  to  the  right- 
hand  side  of  Eq.  (22),  since  integrals  over  z  of  rapidly 
oscillating  terms  can  be  neglected  with  respect  to  in¬ 
tegrals  of  slowly  varying  functions  (rotating-wave  ap¬ 
proximation).  1  More  precisely,  in  order  to  benefit  from 
this  approximation,  one  has  to  assume  that  the  c„’ s  do 
not  undergo  relevant  variations  over  distances  of  the 
kind  [&,(<o)  -  *(»)]*  and  [8,(w)  -  fcU)  -  U')  -  6.(ui)]*\ 

unless,  respectively,  m  «  *  and  m*j,  r*  s  or  m  *  s, 
j*r.  Under  the  same  hypothesis,  it  follows  that  the 
c„’s  do  not  vary  over  toe  distance  [&,(&>)  -  6/ui)  ♦  8,(u/) 

-  6«(w')rl>  unless  m  *  j,  r»  »  or  m  ■  x.  t  =  k  since 

A»(<*>)  -  8|(w)  +  6,(01')  -  6,(ui')“  6»(o>)  -  6<(ui)  +  8,(01)  -  8,(o>) , 

(24) 


f  (a;,(z')a;(0)>  exp {-,  [8,(0.')- 6/01')]  z'}dz',  (22) 

Jo 


provided  that  the  condition  I  o>  -  01' I  /01 « 1  is  fulfilled 
[see  Eq.  (3)].  Equation  (24)  can  be  justified  by  observ¬ 
ing  that 

6,(01)  -  8/(o>)  +  8,(w')  -  6,(01') 

“  6»(o>)  -  6/(oi)+  6,(01)  -  6,(w)  +  (^Bz.  -  (oi'  -  01) , 

(25) 

and  that 

|  (rf&/doi  -  d&, /du>){u>'  -  01)  | 

«  1 6.(01)  -  61(01)  +  ft(oi)  -  6«(o»)l 

Eq.  (26)  being  always  valid  if  toe  differences  between 
toe  inverses  of  the  group  velocities  are  of  toe  same  or¬ 
der  of  the  differences  between  toe  Inverses  of  toe  phase 
velocities. 

By  taking  advantage  of  the  rotating-wave  approxima¬ 
tion,  Eq.  (22)  yields 
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■fcXmrb-  «.  <■»')  =  *„.(*,  ID,  «') 

+  E  CiK*(z')K*rm  exp{-  ifofro')  -  */«/)]«'}*' 

4  •'O 

+  J~  (Km(2'!«»)  dz'),  w  *  r  .  (27) 

It  Is  not  diHicuit  to  show  (see  Appendix)  that  the  co¬ 
efficient  of  X.„  on  the  right-hand  side  of  Eq.  (27)  has  a 
nonpositive  real  part,  so  that  Eq.  (27)  can  be  rewritten 

-j^X.riz,  'o,  u)=-gmrX„(z,  id,  id')  ,  m*r,  (28) 


be  rewritten,  with  the  help  of  Eq.  (11), 

-  ~  £  iP.U,  0)  ,  (33) 

so  that  finally  one  has 
~  iP.U,  /)>  +  j~t  &•>(*•  0) 

=  Z*.A<PA2,  <»-(P.U,  i)>],  (34) 

4 

which  constitutes  the  set  of  coupled  equations  describing 
the  propagation  of  a  pulse  in  a  multimode  lightguide.  At 
this  stage,  distributed  losses  can  be  easily  incorporated 
by  adding  the  term  -  2  a .{P.iz,  /))  on  the  right-hand  side 
of  Eq.  (34). 


with  Regnr^O,  whose  solution  is 

X.r(z,  id,  id')  =  (c.iO,  <d) c? (0,  <d')> „ e’t*r‘  ,  m*r , 

(29) 

where  c„(0,  cd)*c„(z  =0,  id). 

Equation  (29)  shows  that  the  cross-correlation  terms 
in  the  mode-field  amplitudes  practically  vanish  after  a 
distance  of  the  order  1/ R eg., ;  or  that  they  remain  iden¬ 
tically  zero  if  they  possess  this  value  at  z  =  0  (excitation 
by  a  spatially  incoherent  source). 

For  the  diagonal  terms  m  *  r,  we  obtain  from  Eq.  (22), 
<D,  v')  =  “(ZJ  h-lJXmm(z,  <D,  </) 


+  E^-i63^  [&.(<*>)-&.(«')]  Z 
+  i[P,{w)-0j{w))z}Xjj{z,  (D,  w'). 


where 


(30) 


K,=  (K^Ol) 


xexp{i[^J((Do)-ft,((Do)]z'}<fz'  (31) 

is  a  real  non-negative  quantity  (see  Appendix).  In  de¬ 
riving  Eq.  (30)  we  have  taken  advantage  of  the  fact  that 
the  value  of  h.,  is  not  substantially  modified  by  the  sub¬ 
stitution  id  -  w0  and  id'  -  <d0  performed  in  writing  Eq. 
(31). 


The  set  of  Eqs.  (30)  describes  in  full  generality  the 
evolution  of  the  autocorrelation  functions  of  the  mode- 
field  amplitudes  and  allows  in  principle  the  evaluation 
of  <Pt(z,  /))  by  means  of  Eq.  (14).  It  is,  however,  more 
convenient  to  deal  directly  with  a  set  of  equations  con¬ 
necting  the  average  powers  (P.iz,  ())  carried  by  each 
mode.  By  recalling  the  expression  oi  .»’„iz,  ()>  given  by 
Eq.  (14),  one  obtains,  with  the  help  of  Eq.  (30), 


£  iP.U,  t))  =  -  (g  h.]j  (P.iz,  /)) 

♦  E  *.4<^4<*.  <»  *P  f  du»  f  dw'  i[8,(id')  -  8.(id)] 

|#W  J 

xexpl<(w-w')<]exp{-ilA.(w)-d.(w')]z}X,J(z,  id'). 

(32) 

The  last  term  on  the  right-hand  side  of  Eq.  (32)  can 


The  set  of  Eqs.  (34),  which  describes  the  time-depen¬ 
dent  propagation  of  the  average  modal  powers,  is  iden¬ 
tical  with  the  one  worked  out  by  previous  authors. 1-1 
The  present  derivation  has  the  advantage  of  clearly 
showing  that  the  limits  of  validity  of  Eq.  (34)  rely  on  the 
small-bandwidth  assumption  expressed  by  Eq.  (3),  as 
well  as  on  the  hypothesis  of  slow  variation  of  the  field 
amplitudes  required  also  in  the  stationary  case. 1,4 

IV.  PROPAGATION  OF  A  FREQUENCY-MODULATED 
SIGNAL 

The  form  of  the  set  of  Eqs.  (34)  implies  that  the  aver¬ 
age  power  (P.iz,  /))  carried  by  each  mode  is  not  in¬ 
fluenced  by  the  frequency  modulation  of  the  input  signal; 
that  is,  that  this  kind  of  modulation  cannot  be  transmitted 
by  the  fiber,  as  long  as  one  collects  all  the  light  coming 
out  from  the  whole  fiber  section.  This  circumstance  is 
no  longer  true  if  one  considers  the  power  transmitted 
through  a  smaller  area  o  of  the  section.  In  this  case, 
the  orthogonality  of  the  modes  is  no  longer  effective  in 
canceling  the  cross-correlation  terms  between  dif¬ 
ferent  modes,  and  the  expression  of  the  power  (P")  car¬ 
ried  through  this  area  contains  also  nondiagonal  contri¬ 
butions.  More  precisely,  one  has,  according  to  Eq. 

(5), 

(P?)=^Re  ££  T.r,  (35) 

with 

Tmr(z,  0=  JJ  dxdyE.it,  (D0)*  Hj  (r,  (D0)-e. 

0 

*  |  dw  j" du>'  exp {j(id  -  id') l 

-  ttA.((D)-  8,(id')]z}X.,(z,  id,  id')  .  (36) 

Recalling  Eq.  (29),  one  can  rewrite  the  nondiagonal 
contributions  as 

T.r(z,  <)«exp(-g„,z)  jj  dxdyE.it,  iD0)*H*(r,  id0)  •  e, 

x  J  dlD  J  dw'  exp  {»(id  -  o')  I- i  [d«(w)  -  dr((D')]  z} 
X<C«(0,  id)c*(0,  id'))„  ,  m*r,  (37) 
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which  shows  that  the  Tmr’s  ( m*r )  evolve  as  in  the  ab-' 
sence  of  mode  coupling,  except  for  the  fact  that  their 
envelope  undergoes  an  exponential  decay  with  the  fiber 
length  z.  In  particular,  if  the  source  is  spatially  in¬ 
coherent,  then  Tmr(z,  f)  =  0  (m*r),  and  (P*>  is  not  in¬ 
fluenced  by  the  modulation  of  the  signal.  If,  conversely, 
the  source  possesses  some  degree  of  spatial  coherence, 
then  the  Tmr’s  ( m*r )  are  sensibly  different  from  zero 
over  a  certain  fiber  length,  and  Eqs.  (9)—  (12)  can  be 
recalled  for  determining  their  evolution,  thus  getting 

T”C-  ('('-tH-#.’ 

(38) 

If  one  assumes  the  exciting  source  to  be  a  frequency- 
modulated  single-mode  laser,  then 

F (t )  =  exp  [t<>(f)  +  »>/*] ,  (39) 

where  <p(t)  is  a  normally  distributed  random  function 
responsible  for  the  laser  linewidth,  y  is  a  real  constant 
proportional  to  the  frequency- modulation  rate,  and 

T„(z,  t) exp(- gmr  z ) exp  1 1/V»  ~  1'/Vrl 

(40) 

fc=  1/6 v  being  the  coherence  time  of  the  source.  Equa¬ 
tion  (40)  shows  that  Tmr  goes  to  zero  over  the  shortest 
of  the  two  distances  1  /g„  and  tcl  l/Vm-l/VTl,  the  lat¬ 
ter  being  the  traveled  length  over  which  two  modes  m 
and  r  started  from  z  =  0  at  the  same  instant  acquire 


a  time  delay  larger  than  tc .  *  On  the  fiber  length  over 
which  the  Tmr’ s  ( m*r )  are  still  different  from  zero, 
their  time  dependence,  which  is  unaffected  by  mode 
coupling,  determines  the  time  behavior  of  {P\),  since 
the  diagonal  terms  Tmm,  bfcing  proportional  to  (Pm),  do 
not  depend  on  time. 

V.  BEHAVIOR  OF  FLUCTUATIONS  IN  THE 
NONSTATIONARY  CASE 

In  order  to  deal  with  practical  situations,  in  which  a 
single  fiber  is  usually  employed,  it  is  desirable  to  eval¬ 
uate  the  fluctuations  of  the  relevant  quantities  around 
their  average  values.  In  the  monochromatic  case,  this 
has  been  done1'4  by  calculating  for  a  lossless  fiber  the 
asymptotic  value  (in  z)  of  the  normalized  variance 

((F1*)  —  (■P»>i)/(fJ*>i  .  (41) 

For  large  z,  the  resulting  value  is  1,  which  implies 
a  100%  uncertainty  on  the  actual  value  of  the  mth  mode 
power  in  a  single  fiber.  This  rather  disappointing  situ¬ 
ation  drastically  changes  in  the  case  of  a  stationary  sig¬ 
nal  possessing  a  finite  bandwidth  6( o,  where  it  has  been 
shown  that  the  asymptotic  normalized  variance  tends  to 
vanish. '  In  this  section,  we  study  the  evolution  of  the 
fourth-order  field -amplitude  correlation  functions  rela¬ 
tive  to  different  modes  and  different  frequencies,  which 
will  allow  us  to  generalize  the  result  of  Ref.  9  to  non¬ 
stationary  situations. 

The  starting  point  is  the  set  of  Eqs.  (19),  which  im¬ 
mediately  furnishes 


2j(c«U,  <*>)e»(z,  <*>*))  „(c,(z,  «/')c*(z,  u>" ')), 


=  <*>) (fj(z,  <*>)c*(z,  <*)'))„  <c,(z,  w")c* (z,  </"))„  +  A*,(z,  u)')(c„(z,  w)c*(z,  w')>„(c,(z,  u>")c*(z,  u/")>» 


+  ^(z,  w")(c,(z,  oi)c*(z,  w'))„(c/(z,  w" )c* (z,  oi'"))„  + A*/z,  u>")(cM(z,  u >)cj(z,  u")cj (z,  a/")).}. 

(42) 

It  is  worthwhile  to  emphasize  at  this  point  that  the  source-average  (  )„  in  Eq.  (42)  operates  on  mode-amplitude 
pairs  successively.  In  fact,  the  basic  quantity  measured  by  the  detector  in  a  single  fiber  is  the  modal  power  PJz,  t) 
[see  Eq.  (8)],  which  contains  the  second-order  average  (c„(z,  w)c*(z,  <»')>„,  so  that  the  evolution  of  ensemble  aver¬ 
ages  of  the  kind  <P„(z,  t)Pn(z,  t ))  involves  the  fourth-order  products  <c„(z,  u>)e*(z,  u >'))„(c„(z,  w')cj(z,  w*)>„.  in 
turn,  the  evolution  of  (c„(z,  <*>)c*(z,  <«'))„  is  connected  with  nondiagonal  second-order  products  of  the  type 
u))c;(z,  (*)'))„,  which  justifies  the  appearance  of  <  )„<  )„  in  Eq.  (42). 

By  taking  r-m,  s  =  n  and  v  =  u>',  =  it  is  not  difficult  to  obtain  from  Eq.  (42),  by  means  of  a  direct  ex¬ 

tension  of  the  procedure  outlines  in  Sec.  m,  the  following  set  of  equations: 

~tP.(z,  w)£„(z,  w")>  =  -^£  hmj^{P„(z,  <ji)Pn(z,  '/')>  +  £  u)Pn(z,  ui")) 

+  (<(<*(*»  w")cm(z,  (*>"))„ (c„(z,  u.>)c*(z,  u))>„>exp{f[d„(w)-^(w")  +  8B(u)")-8„(w)]z}(l-6w„)  j  (K'-„(z')A-n(0)) 

x  exp  {-<[/»>'') -/UOM<f*'  +  c.cW(£  «^(*.  '■>)&*,  w)>.<c«(z,  'u")cj(z,  w")>J 

xexp{f  [A,(w)  -  +  |B>")  -  ff,(w)]  z}6mt  <*„,{*'>*:,«>)>  exp{- i [(*,(»”)  -  &,(«">]  z'}d*'  ♦ c.  c)  *Nm  , 

(43) 
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where  Nn  represents  same  terms  with  the  exchange 
m  —  n,  ui  — w",  having  defined 

£.(*.«)  =  < |c„U,  a>)|*>„.  (44) 

While  the  set  of  Eqs.  (43)  possesses  a  general  validity 
under  the  usual  hypotheses  justyfying  the  rotating-wave 
approximation,  it  can  be  cast  in  a  particularly  simple 
form  if  the  relation 

|(^- S')(w-u’")lz>>2,r’  j*m  (45) 

is  fulfilled  and  if  the  field  amplitudes  do  not  appreciably 
vary  over  a  distance  z‘m,i>  given  by  [compare  with  Eq. 
(40)] 

aj-"  =  |  (1/Vm  -  1/V,)\ *‘  2 ir/Ow  .  (46) 

In  fact,  in  this  case,  one  can  again  take  advantage  of 
the  rotating-wave  approximation  on  the  right-hand  side 
of  Eq.  (43),  which  reduces  to 

^<P.(2,  wAU,  «"» 

=  32^1  [C^iU,  <*>)•?,,(«,  U>"))  -  {Pjz,  «j)P„(z,  «)")>] 

i 

+  «"AU,  «»-(?.(«,  «")P.U,  <*>)>]  • 

(47) 

From  Eq.  (47),  in  the  case  m  =  n  and  for  values  of  z 
large  enough  so  its  left-hand  side  vanishes,  it  is  im¬ 
mediate  to  deduce  the  following  asymptotic  equality: 

tf,U.  «)#.U.  «"».♦  <*,(«,  «». 

=  2(P.U,  w)P.(a,  «"».  ,  (48) 

where  the  subscript  a  stands  for  asymptotic.  On  the 
Other  hand,  energy  conservation  along  the  fiber  requires 
that 

32*M,  *>)=  32  PA0,  w)  for  every  w,  (49) 

j  / 

which  in  turn  implies 

32  «>"» +\32232&m,  *)pm u,  <*>")> 

m  *  m*j 

3232&M,  »")pjz, «» 

w>Z^(°*  <*>").  (so) 

M  i 

a  relation  which  asymptotically  yields,  with  the  help  of 
Eq.  (48), 

£(£„(*,  w)P„(x,  «")>.+Af(JV-l) {PJz,  *)PmU,  </')>. 

s 


=  32  £»(0>  «>  E  ^(0,  u") ,  '  (51) 

■  ) 

where  N  is  the  total  number  of  guided  modes.  Recalling 
now  the  asymptotic  equipartition  condition,  expressing 
the  fact  that,  in  a  lossless  case,  the  power  per  mode  is 
the  same  in  every  mode  for  z  sufficiently  large, 1 

< PmU ,  (*»))«=  jj  32  PA°f  u>)  for  every  u > ,  (52) 

Eq.  (51)  can  be  rewritten 

32<£mU,  «)P.(t,  y'))a+N(N-  l)(Pm(z,  w )PmU,  w"))« 

m 

=  Nt(Pjz,  *))a<Pjz,  »")).,  (53) 

from  which  follows  the  relevant  expression 

(Pm(z,  u)Pm(z,  u>"))a  =  (Pjz,  to >)),(pm(z,  <u")>fl  .  (54) 

We  now  observe  that  the  energy  Im(z)  pertaining  to  the 
mth  mode 

Uz)m£pjz,tUt,  (55) 

can  be  expressed  by  means  of  Eqs.  (8)  and  (44)  as 

/„(*)  =  2tt.P  J  Pjz,  w)du>  .  (56) 

The  normalized  variance  of  the  statistical  quantity  /„ 
can  be  accordingly  written 

°  (1  du>  J 

~  j  dui  j  du>’r{Pm(z,  u))>C P„(z,  w")>) 

x  Q  d*  |  dw"  (Pm(z,  -o)>  <Pm(z,  w")>)  *  .  (57) 

Since,  for  x  »  z{cm’t\  most  couples  of  frequencies  in 
the  relevant  integration  domain  fulfill  Eq.  (45),  Eq.  (54) 
applies  and,  consequently,  the  normalized  variance 
tends  asymptotically  to  vanish,  that  is 

W*>.  -<4)!)/<4>J-0.  (58) 

This  Implies  that  no  statistical  uncertainty  is  present 
in  /„  so  that  equipartition  of  the  total  energy  among 
the  various  modes  takes  place  in  the  single  fiber.  As 
a  particular  case  of  this  general  statement,  the  asymp¬ 
totic  power  Pm  has  the  same  property  in  the  stationary 
regime,  as  it  has  been  shown  elsewhere,  *  since  in  this 
situation  lm  and  Pm  are  proportional. 

The  complete  investigation  of  the  behavior  of  fluctua¬ 
tions  in  the  nonstationary  regime  requires  the  analysis 
of  the  quantities  (Pi(z,  /)>  and  (Pm(z,  l)Pn(z,  /)>.  To  this 
end,  the  system  of  Eqs.  (42)  has  to  be  specified  to  the 
case  r*m,  s*n,  since 


(P.U,  t)P,(t,  t))*P*  j  dw  j  dw'  |  du>"  J  du>"'  exp[i(u>  -  w')f]  exp  {i(w"  -  </")<  -  »[6,,(u>)  -  &.(“')]  z 

-4ft.(u") )]x}«c.U,  w)c-U,  w')).<c.U,  w"kSU,  «'")>■>  •  (59) 

By  means  of  the  general  procedure  followed  in  this  paper,  it  ia  possible  to  deduce  a  set  of  equations  governing 
the  evolution  of  ((c„(s,  w )<£(*,  w')>„\c.(*,  w")c*  (*,  <u' "))»>.  In  order  to  simplify  the  notation,  let  us  define 
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Ymr(t,  ",  u>')  =  exp[i0,(u)')z -i0„(u))z]<c,.(z,  u)c*(z,  u>'))„  . 
The  tinal  form  of  the  closed  system  of  equations  then  reads 


(60) 


'o,  u>')Ym(z,  uj",  a,'"))  =  -  /3„(w)  +  0>'")  -  &„(*")}  (Ymm(z,  u>,  u,')  YJz,  ui",  a,'")) 


^  »')Y„U,  ",  "')*«(*,  u>'"» 

+{Ymn(z, «,  U)')y„.(z,  •„",  hrtL<y.«(Z, ",  "')y«»(2, «',#»+<r-u, ",  «')**.<*, "'")>])  («« 


and 

d 


^<y«U,  w,  </)y,,(z,  <*>",  w'M)>=*U("')-3-(")+^(""')-d,("")](y«rfe,  ",  "')ymU, "*)> 

E  "')yn.U,  "",  u"))  +  T,h«(Ymk{z,  «,  *')Y>m(z,  u",  «'"»  +  2  *^<y,T(z,  w,  «') 

\*#r  w«  /  k#r  Jt#m 


x  y^U,  "”,  -  Kt  [<yrr(z,  ",  0)')y„,,(z,  to",  u’"))  +  (Y„m(z,  U>,  0>')y,,(z,  U>",  «'")>]  ,  m*  T  . 


(62) 


r 


The  set  of  Eqs.  (61)  and  (62)  represents,  together  with 
the  set  of  Eqs.  (30),  a  fairly  complete  analytical  descrip¬ 
tion  of  nonstationary  propagation.  In  fact,  while  Eqs. 

(30)  concern  the  evolution  of  the  average  power  per 
mode,  Eqs.  (61)  and  (62)  allow  one  to  describe  the  fluc¬ 
tuations  aroand  this  quantity  in  terms  of  fourth-order 
averages.  One  can  observe  that,  in  the  general  case,  it 
is  not  possible  to  obtain  a  system  of  equations  containing 
only  the  quantities  (P.P*),  while  this  is  possible  in  the 
monochromatic  stationary  propagation. 4  This  is  rela¬ 
ted  to  the  essential  role  of  the  average  over  the  source, 
which  does  not  allow  the  nondiagonal  terms  (Ymr(w,  to') 

*  y^(w",  to"'))  to  be  rearranged  as  products  {Ymm{u, 

<o'")y,T(<o",  u,')). 

VI.  CONCLUSIONS 

We  have  derived  a  set  of  statistical  coupled  equations 
describing  the  time-dependent  propagation  of  a  finite - 
bandwidth  carrier,  which  can  be  either  amplitude  or  fre¬ 
quency  modulated,  placing  into  evidence  the  limits  of 
validity  of  the  approach,  which  involve  both  the  value  of 
the  bandwidth  and  the  variation  of  the  mode-amplitude 
coefficients.  In  particular,  the  time  evolution  of  a  fre¬ 
quency-modulated  signal  does  not  depend  on  mode  cou¬ 
pling,  thus  possessing  a  behavior  completely  different 
from  that  of  an  amplitude-modulated  signal. 

We  have  set  in  a  rigorous  form  a  theorem  describing 
the  smoothing  of  mode-power  fluctuations  over  the  fiber 
ensemble,  due  to  the  finite  bandwidth,  so  as  to  include 
mode-energy  fluctuations  in  r.onstatlorarv  situations. 

The  general  case  of  time -dependent  propagation  of 
fluctuations  has  been  also  considered,  and  a  closed  sys¬ 
tem  of  equations  enabling  to  describe  it  has  been  de¬ 
rived. 

APPENDIX 

If  one  defines 

/ 1  \i/*  /■*/* 

■F*,J  \l)  (Al) 


one  can  write 

j  rLIZ  rL/l 

Ip’ll*- T  I  *'  j  dz(KMi(z)K*U'» 

J-L/t  J-L/l 

x  exp[f(0,  -  0  J(z  -  z')] ,  (A2) 

so  that,  remembering  Eq.  (31),  and  taking  advantage  of 
the  finite  range  of  the  correlation  function 

(K./zJOz'))  =  <#*,(*  -  z')C(°)>  ,  (A3) 

one  has 

*-,=  lim|Pj*aO.  (A4) 

Let  us  now  consider  Eqs.  (27)  and  (28).  Since  the 
quantity  (Kmt{z)Kjm{0))  is  a  real  even  function  of  z,  one 
can  write 

Re?-,  =  \  Z  A-i  +  i  Yj  hlr 

*  i*m  *  i*r 

*  1  is  <i(r  • 

(A5) 
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The  effect  of  strong  mode  coupling  on  modal  dispersion  in  optical  fibers  has  been  investigated.  The 
pulse  dispersion  turns  out  to  be  qualitatively  different  from  the  one  relative  to  the  weak-coupling 
case,  while  it  exhibits  a  drastic  reduction  as  compared  with  that  of  the  uncoupled  case.  The  role  of  the 
initial  pulse  length  and  of  the  source  coherence  time  has  been  elucidated. 


I.  INTRODUCTION 

The  role  of  mode  coupling  in  the  propagation  of  guided 
modes  in  a  multimode  lightguide  is  well  established.  It  can 
he  negative  or  positive,  according  whether  the  purpose  is  to 
let  the  lower-order  modes  propagate  without  progressively 
sharing  their  energy  with  the  higher-order  ones,  or  to  reduce 
the  modal  dispersion  of  an  optical  fiber.  While  the  first 
statement  can  be  easily  understood,  the  second  one  is  not  so 
obvious  and  has  to  be  proved  analytically.  In  fact,  the  existing 
theory12  is  based  on  a  statistical  approach  and  it  describes  the 
propagation  of  the  average  power  in  the  mth  mode  (Pm(z,t)), 
where  the  averaging  operation  < . )  is  meant  to  be  per¬ 

formed  over  an  ensemble  of  many  macroecopically  equivalent 
lightguides  differing  among  themselves  be.^se  of  micro¬ 
scopically  random  imperfections. 

In  this  frame  it  can  be  shown  that,  undei  suitable  assump¬ 
tions  on  the  coupling,  the  centers  of  mass  of  the  wave  packets 
describing  the  <Pm(z,t))'s  tend  to  travel  with  a  common  av¬ 
erage  velocity,  while  the  pulse  widths  turn  out  to  be  propor¬ 
tional  to  (o  +  bz)l/2,  a  and  b  being  two  prsitlv?  constants,  and 
2  the  traveled  distance.  This  fact  is  interpreted  as  a  mecha¬ 
nism  of  reduction  of  the  modal  dispersion  in  optical  fibers,1-2 


since  one  would  expect,  should  the  modes  travel  indepen¬ 
dently  as  they  do  in  the  absence  of  mode  coupling,  a  distortion 
proportional  to  the  traveled  distance  z. 

The  main  limitation  of  the  statistical  approach  is  connected 
with  the  presence  of  fluctuations  around  the  average  value, 
which,  if  too  ample,  do  not  allow  the  average  result  to  be 
confidently  applied  to  the  single  fiber — with  which  one  most 
often  deals  in  practical  situations.  It  has  been  demonstrated3 
that  in  the  case  in  which  the  fiber  is  excited  by  a  polychromatic 
source  there  is  a  fiber  length  L  <*  1/TC  (Tc  being  the  coherence 
time  of  the  source)  such  that  the  fluctuations  of  the  energy  per 
mode  Im(z,t): 

Im(z,t)  -  J*_"  Pm(z,t)  dt  (1) 

tend  to  vanish  for  z  larger  than  L .  However,  this  result  does 
not  concern  the  behavior  of  Pm(z,t)  itself,  but  for  the  sta¬ 
tionary  case  in  which  Pm  (z,t)  is  time  independent.  Besides, 
the  statistical  approach  does  not  give  a  correct  result  for  the 
single  fiber  when  one  has  to  evaluate  the  cross-correlation 
term  between  the  field  em(r,z,t)  and  e„(r,2,t)  pertaining  to 
the  mth  and  nth  mode,  that  is 
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where  the  symbol  ( . ).v  indicates  the  averaging  operation 

over  the  fluctuations  of  the  source  exciting  the  lightguide.  As 
a  matter  of  fact,  it  has  been  shown*  that  (Tnm )  goes  to  zero 
over  a  distance  not  larger  than  Te  1 1/V„  -  1/ V„  |  ~ 1  ( Vn  and 
Vm  being  the  nth  and  mth  mode  group-velocity),  which  is  just 
the  distance  over  which  it  would  go  to  zero  in  the  absence  of 
mode  coupling.4  This  obviously  cannot  be  the  case  for  the 
single  fiber,  since  the  presence  of  a  coupling  certainly  gives  rise 
to  a  correlation  between  the  modes  over  a  longer  distance. 

The  previous  considerations  make  clear  that  it  would  be 
desirable  to  be  able  to  find,  at  least  for  some  simple  workable 
models,  an  analytic  solution  for  the  problem  of  propagation 
in  a  deterministic  case  (that  is  without  resorting  to  the  sta¬ 
tistical  approach),  in  order  to  compare  it  with  the  above 
mentioned  results.  This  would  also  allow  us  to  evaluate  in 
a  correct  way  the  Tnm 's,  whose  behavior  may  furnish  a  simple 
way  of  gaining  information  on  dispersion — and  thus  on  mode 
coupling*  In  this  paper,  this  program  is  carried  out  intro¬ 
ducing  the  hypothesis  of  strong  coupling,  a  case  interesting 
per  se  which  cannot  be  investigated  by  means  of  the  statistical 
treatment  which  covers  weak  coupling,  and  considering  a 
simple  model  in  which  only  two  modes  interact 

In  the  framework  of  our  approach,  strong  coupling  means 
that  the  magnitude  of  the  coupling  constant  K12is  much  larger 
than  the  difference  |di  -  d2|  between  the  propagation  con¬ 
stants  of  the  two  unperturbed  modes  (see  Gq.  (24)].  Under 
this  assumption,  we  consider  the  situation  of  resonant  cou¬ 
pling,  for  which  the  characteristic  spatial  periodicity  lofKu 
if  of  the  order  of  |di  —  d2| -1,  end  that  of  slowly  varying  cou¬ 
pling,  for  which  l  »  |di  -  /32| ~l-  In  both  cases,  propagation 
is  significantly  affected  by  the  fact  that  |Ki2|  is  large  (while, 
as  is  well  known,  the  fulfilment  of  the  resonant  condition  is 
essential  in  the  weak-coupling  regime2).  Dispersion  turns  out 
to  be  drastically  reduced  with  respect  to  the  case  in  which 
coupling  is  absent,  and  its  qualitative  behavior  differs  from 
what  one  would  expect  according  to  the  statitistical  method, 
in  that  the  pulse  spreading  turns  out  to  be  proportional  to  the 
fiber  length  z  instead  that  to  zin. 

The  analysis  of  the  propagation  in  the  presence  of  strong 
coupling,  besides  allowing  us  to  complete  the  description  of 
the  effects  of  mode  coupling  on  modal  dispersion,  can  be  rel¬ 
evant  for  the  study  of  propagation  in  mode  scramblers,  which 
are  strongly-coupling  fiber  samples  able  to  excite  all  guided 
modes  in  a  repeatable  manner.6 

II.  DESCRIPTION  OF  PROPAGATION  IN 
UGHTGUIDES 

The  transverse  electromagnetic  field  propagating  in  a  cy¬ 
lindrical  guiding  structure  can  be  expressed  as  the  superpo¬ 
sition  of  the  fields  pertaining  to  »8<  h  guided  mode  in  the 
form3,7 

E(r,z,t)  -  Z  em(r,z,t), 

m 

with 

em{rj,t)  -  Em(r)om(r,t),  (4) 

where  the  En(r)'s  represent  the  modes  of  the  ideal  guiding 


structure  (that  is  without  mode  coupling)  and  the  expansion 
coefficients  am(z,t)  are  defined  through  the  relation 

am(*.0  “  J"  cm(2,w)e‘-'-**»<">*dto).  (5) 

In  Eq.  (5),  dm(<*>)  is  the  propagation  constant  of  the  mth 
mode  and  the  cm 's  depend  on  2  because  of  the  presence  of 
mode  coupling  [otherwise  one  would  have  cm(2,b>)  ■  cm(0.u)). 

It  is  possible  to  show6  that  as  long  as  only  forward  traveling 
modes  are  considered  the  transverse  part  of  the  magnetic  field 
obeys  equations  identical  to  Eqs.  (3)  and  (4),  provided  the 
substitution  Em  -»  Hm  is  made.  From  the  above  consider¬ 
ations  it  follows  that,  in  order  to  evaluate  second-order  aver¬ 
ages  of  the  kind 

G(ri,r2,ti,t2,2)  *  <E(ri,2,t1)-E,(r2,2,t2)>.v,  (6) 

or 

P"(z,t)  =  ~Re  f  dr  (E(r,2,t)XH»(r,2,t)).v-el  (7) 

(e,  being  a  unit  vector  in  the  z  direction),  representing,  re¬ 
spectively,  the  mutual  coherence  function  and  the  power 
carried  by  the  electromagnetic  field  through  a  given  area  a  of 
a  fiber  section,  it  is  sufficient  to  investigate  the  behavior  of  the 
quantities 

(om(z,ti)o;(z,t2)  >.v  (8) 

or,  equivalently,  of 

<cm(z,w)c‘(z,w'))„.  (9) 

According  to  Eqs.  (5),  (4),  and  (7)  and  the  mode  orthogon¬ 
ality,  it  turns  out3  that  (am(z,t)a'„(z,t))tv  is  proportional  to 
the  power  Pm  (z,t )  carried  by  the  mth  mode  through  the  whole 
fiber  section,  and  that 

P-(*,t)-ZPm(z,t),  (10) 

m 

while  the  nondiagonal  terms  <am(2,t)a’(z,t))«v  can  be  con¬ 
nected  with  the  degree  of  correlation  between  the  various 
propagating  modes.  In  particular  one  has 

P’U,t)m  J ReZLPmn  (“m (z,t)a’n(t,t)),v 

a  m  n 

*  Re  L  (F’mJF~m)Pm(z,t) 

m 

+  ;ReE  E  F’mn  <,am(z,t)a\(z,t))„,  (11) 

with 

F’mn  «  J  dr  Em  (r)  X  H*„(r)  •  e„  (12) 

showing  that  the  nondiagonal  terms  represent  interference 
contributions  between  the  various  modes,  which  are  always 
present  whenever  a  is  finite  (otherwise  they  disappear  due  to 
the  orthogonality  between  different  modes).  Thus,  the 
analysis  of  P'(z,t),  as  compared  with  P“(2,t),  can  furnish  a 
way  of  estimating  the  degree  of  correlation  between  the 
propagating  modes. 

In  the  absence  of  mode  coupling,  one  has  cm (2, «)  «  cm(0,n>), 
which,  together  with  the  approximate  relation 

dm  («)  “  dm  ("0)  +  (<*'-  O>o)/Vm  (13) 
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(which  implies,  as  a  necessary  condition,  the  ratio  between  the 
bandwidth  of  the  propagating  signal  and  the  central  fre¬ 
quency  <jo  to  be  much  less  than  one),  furnishes  (see  Eq.  (5)] 

am(z,t)  -  om(0 ,t  -  z/Vm),  (14) 

where  Vm  =  (d0m/du)Zl«1  is  the  group  velocity  of  the  mth 
mode.  This  amounts  to  saying  that  Pm  (z,t )  propagates  with 
this  velocity  and  that  the  nondiagonal  terms 

(Sn  (z,t  )fl  „  (z,t )  )BV 

-  (am (0,t  -  z/Vm)a'Mt  -  z/V„))„  (15) 

vanish  whenever  the  fiber  length  is  such  that  they  acquire  a 
time  delay  larger  than  the  coherence  time  Tc  =  It  thus  of  the 
source;  that  is, 

\z/Vm-z/V„\>Tc.  (16) 


III.  PROPAGATION  IN  THE  PRESENCE  OF  MODE 
COUPLING:  NONRESONANT  CASE 


In  order  to  have  an  analytically  solvable  model,  We  consider 
a  simple  case  in  which  only  two  modes  1  and  2  are  interacting. 
The  evolution  with  z  of  the  mode  amplitudes  ci(z,w)  and 
c2(z,u)  is  described  by  the  following  set  of  equations8: 

^^-/f12(z)e^«->*c2(z,w), 

„  i  >  <17> 

-X;2(z)e-^«->*c,(z,a.). 

where  A(oj)  =  0,(u)  -02(ai),/fi2(s)  is  a  z -dependent  coupling 
coefficient  and  the  self-coupling  coefficients  Afu(2).  X22 (*), 
which  would  not  give  rise  to  any  relevant  effect,  have  been 
omitted  for  the  sake  of  simplicity.  By  means  of  the  change 
of  dependent  variables 


ct (z,u>)  *  bi(z,u>)eiMuU/2, 
c2(z,u>)  *  b2(z,u>)e~iMuU/2, 
one  gets  from  Eq.  (17), 

dz 

^-UA/2)b  2--k;25,, 


(18) 


(19) 


from  which,  by  eliminating  b^  one  arrives  at 

b\  -  (Kn/Kn)b\  +  (|X12|2  +  A2/4 

“  i(/Ci2/2X i2)A]  6]  *  0,  (20) 


where  the  prime  indicates  differentiation  with  respect  to  z. 
By  performing  the  substitution 

6i-6iexp(i  j*'(KI2/KI2)dz),  (21) 

one  finally  obtains 

6\  +  IlK.zl2  +  A*/4  -  i(K\i/2KvdA 

+  (K'l3/2Kt2)'  -  </fij/2K12)s)6,  -  0.  (22) 

Let  us  now  assume  Xj2(z)  to  be  a  3iowly  varying  function 
of  z.  This  means  that  there  is  no  Fourier  component  of  the 
coupling  process  that  would  provide  coupling  in  the  sense  of 
first-order  perturbation  theory.  However,  this  fact  does  not 


preclude  an  effective  coupling,  with  a  consequent  influence 
on  dispersion,  for  |Xi2|  large  enough,  as  the  following  deri¬ 
vation  will  make  clear.  The  last  three  terms  inside  the  square 
bracket  in  Eq.  (22)  can  be  neglected,  as  compared  with  |Xi2|2 
+  A2/4,  provided  that  the  condition  l  »  |  1/A|  is  fulfilled,  / 
being  a  characteristic  oscillation  length  of  K\2(z).  The  re¬ 
sulting  equation  can  be  easily  solved  in  the  WKB  approxi¬ 
mation,  under  the  usual  adiabatic  hypothesis 


-7t<|K«l2+*2A01/2l 

az 

|K12|2+  A2/4 


In  particular,  in  the  strong-coupling  regime, 


(23) 


|K,2|2»  a2 


(24) 


and  by  assuming,  for  the  sake  of  simplicity  and  without  loss 
of  validity  of  the  main  conclusions  of  our  derivation,  X12  to 
be  independent  from  Z,  one  obtains 

at(z,t)  *  f  du  exp(—  t(di(w)  +  02(“)]  z/2  +  tut) 

X  (A  exp(t|X12|z  +  iA2z/8|Xi2|) 

+  B  exp(-i |X12|z  -  iA2z/8|X12|)),  (25) 

a2(z,t )  *  i  exp ( — i<p)  J  du  expl— i[8i(ai)  +  /J2(w)]z/2  +  iwt) 

X  (A  exp(i|Xi2|z  +  iA2z/8|X12|) 

-  B  exp(-i|Xi2|z  -  iA2z/8|Xi2|)],  (26) 

with 


A  -  (%)[c,(0.a,)  -  ie>*c2(0,u>)l,  (27) 

B  =  (V2)[c,(0,w)  +  ie‘*c2(0, <*>)], 

and  exp(i»  *  Xi2/|Xt2|.  Assuming  ci(0,a>)  and  c2(0,u>)  to 
possess  the  same  u  dependence,  the  problem  is  reduced  to  the 
investigation  of  the  quantities 

</«(*, 0  *  S  du  ci(0,w)  exp|-  i\0iM  +  /}2(w))z/2| 

X  exp[iut  +  i|Xi2|z  +  iA2(a>)z/8|K12(],  (28) 

and 

Jb(z.t)  -  f  dw  C](0,a>)  exp|— i(di(<*>)  +  /J2(w)(z/2| 

X  expjiwt  —  i|Kj2|z  —  iA2(u>)z/8|Ki2|l,  (29) 

which  is  obtained  from  J\  by  means  of  the  substitution  |KI2| 

--|X12|. 

By  taking  advantage  of  Eq.  (13),  one  can  write 

A(<i>)  *  A(u>o)  +  (o>  —  u>o)/V(-),  (30) 

where 


l/V<-»  -  1/V,  -  1/V2,  (31) 

which,  once  introduced  in  Eq.  (29),  yields 

Jb{z,t)  -  e‘*S  du>  Ci(0,<o)e~‘",Tl  e*w<l  -  z/0b)  (32) 


with 


|diW  +  fe(wo)l*/2  +  wo*/2V<+> -  |K,2|z 
-  A2(wo)*/8|X12|  +  A(wo)wo*/4V«-»|X,2| 

-  <4r/8V<-»,|X12|  (33) 

and 

r2  -  z/8V«-»2|X12|,  (34) 
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having  introduced  the  new  velocities  V{*)  and  0*  defined  by 
the  relations 

l/V™  -  l/Vi  +  1/V2,  (35) 

1/Ob  -  l/2V<+>  +  A(<Jo)/4V<->|/fi2l  -  «o/4W->*|ff,2|.  (36) 

An  analogous  expression  holds  for  </„(z,t)  once  the  substi¬ 
tution  |Kj2|  -»  -|K12|  is  performed,  which  in  particular  de¬ 
fines  another  velocity 

l/0a  -  1/2 V<+)  -  A(«o)/4V<->|/f,2|  +  uJW~*\K12\. 

(37) 

We  have  now  to  consider  the  quantities  JaJ\,  JbJb.  and 
JaJ’tn  and  to  average  them  over  the  fluctuations  of  the  source 
(which,  for  a  stationary  source,  is  equivalent  to  taking  the  time 
average  over  an  interval  of  the  order  of  the  coherence  time  Tc). 
In  fact,  their  linear  combinations  furnish  the  significant 
quantities  (|ai(z,t)|2),v.  (|o2(2,t)|2),vand  <ai(z,t)aj(z,t)),v. 
Thus,  according  to  Eq.  (32),  one  has  to  evaluate  the  quanti¬ 
ty 

<ci(0,w)c;t0,«')>.v  (38) 

This  can  be  done  by  assuming  the  electric  Held  at  z  =  0  to 
possess  a  temporal  behavior  of  the  kind 

e‘“® <F(t)S(t),  (39) 


where  F(t)  is  a  rapidly  varying  function  accounting  for  the 
source  bandwidth  while  S(t)  represents  the  (usually)  much 
slower  amplitude  modulation  of  the  carrier.  By  assuming 

S(t)  .  expi-tVT*)  (40)  J 

and 

«  exp  I-(t'  -  t')2/T?  1.  (41) 

where  the  subscripts  p  and  c  stand,  respectively,  for  pulse  and 
coherence,  one  has  (see  Appendix) 

<ci  (0,u)cf(0,n>')),v 

«  exp(-ft2T?/4  -  S2'2r?/4  +  Qft'7l/2),  (42) 

where 

71  «  r2(T2  +  TAHITI  +  7?)  (43) 

and 

71  -  T*/(2Tl  +  T2),  (44) 

with  ft  =  id  —  uq  and  ft'  *  w'  —  wo-  In  particular,  the  mono¬ 
chromatic  and  the  stationary  cases  are  respectively  recovered 
by  letting  Tc  and  Tp  become  infinitely  large. 

By  taking  advantage  of  Eq.  (42),  one  can  write 


S  do,  S  du,'  exp  [— «2(«t2  +  Tf/4)  -  w'2(-ir2  +  71/4)) 

X  exp  |w[if*  +  w0(T?  -  T\)/ 2)  +  w'(-it*  +  wo (Tf  -  T22)/2  +  w7l/2)|  exp[-w$(7j  -  7l)/2),  (45) 


with  tbm  t  —  z/0b-  After  performing  the  two  integrals  over 
u  and  w',  Eq.  (45)  furnishes 


<«V*).v  “  (l/T)exp[-(t  -z/i’*)2/T2] 

(46) 

where 

T-  (71/2  +  8r V7D1'2, 

(47) 

l/ub  «  l/2V,+)  +  A(«o)/4V<-»|/f12|, 

(48) 

with 

71  -  T\  -  71  -  T2T2/(2T2  +  T2), 

(49) 

Tl  -  71  +  Tf  -  (2Tp  +  T\T])K2T\  +  T\). 

(50) 

An. analogous  expression  holds  true  for  <«/a«/j)tv, 

reads 

which 

«  (1  IT)  exp(-(t  -  z/u.)2/7’2). 

(51) 

where 

l/o.  -  l/2V<+>  -  A(wo)/4V-'|K,2|. 

(52) 

One  can  observe  that  the  group-velocities  va  and  u*  could 
have  been  also  directly  obtained  by  taking,  respectively,  the 
derivative  with  respect  to  u  (evr'uVed  it  «  *  c^)  of  the 
propagation  constants  of  the  fust  and  second  terms  contrib¬ 
uting  to  Oi  (Eq.  (25)]  or  a2  (Eq.  (26)}. 

The  term  <  •/.«/*)„  cannot  be  expressed  in  the  form  of  a 
wave  packet  propagating  with  a  definite  velocity,  and  will  be 

discussed  later  on. 


The  expression  of  (J„«0.v  and  ( </**/*  >»v,  as  given  by  Eqs. 

(46)  and  (51),  shows  that  the  signal  inside  the  fiber  evolves  in 
two  distinct  pulses,  traveling  with  two  slightly  different  ve¬ 
locities  ua  and  vb.  Besides,  the  temporal  width  of  each  pulse 
increases  with  z,  so  that  one  has  to  take  into  account  two  dif¬ 
ferent  sources  of  dispersion.  In  practice,  this  second  kind  of 
dispersion  can  be  neglected  for  the  values  of  z  such  that  (see 
Eq.  (47)] 

T2(z)/T3  <  7V4  a  Tp/4.  (53) 

The  other  contribution  can  be  put  in  a  quantitative  form 
by  introducing  the  time  delay  Td  between  the  centers  of  mass 
of  the  two  pulses  at  a  distance  z;  that  is, 

rrf(z)  -  zd/o*  -  l/oa)  -  A(«„W2V<->|Kl2|.  (54) 

It  is  immediate  to  see  that  this  kind  of  dispersion  is  domi¬ 
nant — the  distance  between  the  centers  of  mass  of  the  two 
pulses  increasing  more  rapidly  than  their  widths.  In  fact,  it 
follows  from  Eqs.  (34)  and  (54)  that 

T4TJ r2  -  4V<-»A(w0)7’3 

■  4u)07’3V(->A(wo)/wo  “  4u»0T3,  (55) 

having  taken  into  account  the  circumstance  that  the  difference 
between  the  inverse  of  group  and  phase  velocities,  V<->  ancT 
wo/  A(wo),  are  of  the  same  order*  Observing  that  Ts  is  of  the  ✓ 

order  of  the  coherence  time,  the  ratio  between  Td  and  r*/T3  Vi 
turns  out  to  be  wq/4w  »  1. 
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Returning  now  to  the  problem  of  evaluating  the  influence 
of  the  term  it  may  be  noted  that  its  importance 

tends  to  become  negligible  for  distances  such  that 

Td(z)  >  Tp,  (56) 

after  which  the  two  pulses  do  not  overlap  anymore. 


IV.  PROPAGATION  IN  THE  PRESENCE  OF 
MODE  COUPUNG:  RESONANT  CASE 

The  resonant  case  corresponds  to  a  sinusoidally  modulated 
coupling  constant 


Ki2(z)  “  2 K  cos(xr). 


(57) 


with 


X  *  A(«o),  (58) 

the  (resonance)  frequency  u>o  having  been  assumed  to  coincide 
with  the  central  frequency  of  the  exciting  source.  By  intro¬ 
ducing  Gq.  (57)  into  Eqs.  (17)  one  obtains 

dc  1(2,01) 


dz 

dc2(z,o>) 

.  3 

dz 


,  Ke-ix,+'^<“*zc2(2,u), 


— K*e‘^-iMu),Ci(z,u), 


(59) 


having  neglected  in  the  right-hand  side  of  Eqs.  (59)  the  terms 
containing  the  rapidly  oscillating  factors  exp|±x'(A(wo)  +  x]*l 
with  respect  to  the  slowly  varying  terms  containing 
exp|±i[A(uo)  -  x]*l.  whenever  |K|  «  |x  +  A(o>)|.  By  using 
this  approximation,  it  is  possible  to  write 

Kl2-  Ke-'*'*'*.  (60) 

One  can  now  use  the  procedure  of  Sec.  Ill,  thus  obtaining 

Oi(2,t)  »  exp  ^  »A(«o)»^  j 


.  fllM  +  02(<n) 
2 

iA2(o>)z 


X  exp^-i 
<  exp  ^i|IC| 


2  +  iuit'j 

_ tA(o)o)A(ai)2 

•1*1  4|K| 

tA2(toi)2 


8|/f| 


.  iA((i>o)A(o))z  iA2(uo)2\  _ 

+— W\  (  ) 


a2(z,t)  **  i  exp 


>A(u)q)2^ 


/  dw 


(*“  2 


2  + 


iwt^j 


where exp(i*)  ■  K/\K\. 

The  expressions  of  a, (2,  t)  and  a2(z,t)  show  that  the  main 
properties  of  propagation  can  be  deduced  following  the  deri¬ 
vation  of  Sec.  III.  In  the  present  resonant  case,  the  group- 
velocities  va  and  t>6  of  both  contributions  to  aj  (and  a2) 
coincide,  since 

1  m  _d_  (d\(o>)  +  fl2(n)  ^  A2(a>)  A(h»o)A(ii))\ 

Va,kmdo)\  2  T  8|K|  *  4|K|  / w*wo 


so  that  the  signal  appears  in  the  form  of  a  single  pulse 
undergoing  a  temporal  broadening  T  given  by  Eq.  (47),  which 
for  large  2  reduces  to 

r*81/2  — *  — (65) 

In  other  words,  the  two  distinct  pulses  of  the  nonresonant 
situation  overlap,  so  that  the  only  cause  of  dispersion  is  given 
by  their  common  broadening. 

V.  INFLUENCE  OF  MODE  COUPLING  ON 
DISPERSION 

One  can  compare  the  results  of  the  previous  sections  with 
the  one  obtained  in  the  frame  of  the  statistical-coupled  power 
theory,  according  to  which  the  temporal  width  of  the  pulse  T, 
is  proportional  to  the  square  root  of  z,  that  is10 

T.  -  (2zF2/V<-»h>'2,  (66) 

where 

h  -  J  d{  <K M)K ;2(0)  >e‘A<"o,f  (67) 

the  angular  brackets  indicating  the  averaging  operation  over 
an  ensemble  of  many  macroscopically  similar  fibers.  Equa¬ 
tions  (54),  (65),  and  (66)  show  that  the  results  of  the  weak-  and 
strong-coupling  cases  are  qualitatively  different  In  partic¬ 
ular,  the  deterministic  approach  brings  under  consideration 
the  role  of  the  coherence  time  Tc  of  the  exciting  source  and 
of  the  pulse  duration  Tp.  It  provides  a  dispersion  which  ex¬ 
hibits  the  same  linear  z  dependence  of  the  uncoupled  case,  for 
which  the  width  Tu  of  the  signal  carried  by  the  two  uncoupled 
modes  is  given  by10 

Tu  -  z/V<‘>,  (68) 

but  which  is  quantitatively  very  different.  In  fact,  one  has, 
Td/Tu  -  A(«o)/2|Ki2|  «  1  (69) 

and 


A  .  .....  A-..,.  .  >A2(«)2  i'A(ujo)A(u>)z  .  iA2(u>o)z\ 

'  ^  «,  +  -JjJj - 

„  /  iAHw)z  .  «A(«o)A(w)z  iA2(«o)*\  1 

Bexp(-,|K|,— ^4-— ^ 


A  -  <%)[ci<0,»)  -  ie'*c2(0,«)J,  (63) 

B  -  (Vi)(c,(0^)  4-  ie'*c2((Ml, 


T/Tu  «  TJTU  -  A(«o)/2|/fl2|  «  1,  (70) 

where  the  first  inequality  in  Eq.  (70)  follows  from  Eq.  (55),  in 
agreement  with  the  intuitive  statement  that  the  presence  of 
a  strong  coupling  must  considerably  reduce  dispersion,  in  both 
resonant  and  nonresonant  situations. 

We  wish  to  stress  again  the  fact  that  the  present  deter¬ 
ministic  method  deals  with  a  case  in  which  the  cm’ s  consid¬ 
erably  vary,  in  the  nonresonant  case,  over  a  distance  of  the 
order  l/A(«o)  (strong  coupling),  while  the  opposite  hypothesis 
is  introduced  in  order  to  deal  with  the  statistical  approach.3 
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Actually,  the  strong -coupling  hypothesis  has  as  a  consequence 
that  the  powers  per  mode  Pt(z,t)  and  Pfa.t)  possess  the  same 
temporal  evolution.  Accordingly,  two  initially  correlated 
modes  travel  without  acquiring  any  mutual  time  delay,  so  that 
(ai(z,t)ai(?,t))|V  does  not  vanish  after  the  distance  defined 
by  Eq.  (16),  which  implies  that  the  transverse  spatial  corre¬ 
lation  between  them  is  preserved  over  a  very  long  traveled 
path. 


VI.  CONCLUSIONS 

We  have  treated  the  propagation  of  a  pulse  of  given  initial 
duration,  injected  into  the  fiber  by  a  source  with  definite 
temporal  coherence  properties,  by  adopting  a  simple  model 
in  which  only  two  modes  are  considered.  This  has  been  per¬ 
formed  under  the  assumption  of  strong  coupling,  by  means 
of  a  deterministic  approach,  which  is  able  to  cover  both  the 
resonant  and  far-from-resonance  cases. 

Strong  coupling  is  effective  in  both  resonant  and  nonreso¬ 
nant  conditions  and  affects  qualitatively  and  quantitatively 
the  behavior  of  dispersion.  Far  from  resonance,  dispersion 
is  associated  with  the  breaking  up  of  the  initial  signal  into  two 
distinct  pulses  proceeding  with  different  velocities,  such  that 
the  corresponding  mutual  delay  exceeds  the  temporal 
broadening  of  the  single  pulse.  Conversely,  this  broadening 
turns  out  to  be  the  only  sburce  of  dispersion  in  the  resonant 
case,  where  the  signal  travels  as  a  single  pulse. 

In  both  cases,  dispersion  is  drastically  reduced  with  respect 
to  the  uncoupled  situation,  and  the  main  qualitative  difference 
with  respect  to  the  weak -coupling  regime  consists  in  the  linear 
dependence  of  pulse  dispersion  on  the  traveled  length. 
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APPENDIX 

According  to  Eqs.  (5)  and  (39),  one  has 

ci(0,<i>)  «  dt'  e--'  S(t')F(t')  e*"*1',  (Al) 


so  that 

(C)(0,w)cj(0,<i>'))*v  «r  J*“  dt'  J’  dr  S(t')S*(t') 

X  e~‘t“'+‘u','<,F(t')F*t'))„,  (A2) 

where  fl  “  a)  -  u»o  and  O'”  <tf'-aio.  After  introducing  the  time 
Fourier-transform  S(f)  of  the  slowly-varying  amplitude  Sit), 
one  can  take  advantage  of  the  stationariety  of  the  rapidly- 
varying  part  F(t) 

<F(t')F*(t')).v  *  G(t'  -  tr),  (A3) 

thus  being  able  to  rewrite  Eq.  (A2)  in  the  form 

(ci(0,<D)et(0,u)')>«v  *  J*  df 

X  J*  *  dr  S(f)S(f  +  «'  -  n)G(T)e-'<u-'»'.  (A4) 

By  recalling  that,  according  to  Eqs.  (40)  and  (41),  one 
has 

S(/)  >  exp(-/2Tp/4)  (A5) 

and 

J"  dr  exp[— i(B  -  f)r\  G(t)  a  exp[-(Q  -  f)2T*/ 4], 

(A6) 

it  is  possible  to  perform  the  integration  in  Eq.  (A4),  thus  ob¬ 
taining  Eq.  (42). 

*On  leave  of  absence  from  Fondazione  Ugo  Bordoni,  Roma,  Italy. 
*S.  D.  Personick,  “Time  Dispersion  in  Dielectric  Waveguides,”  Bell 
Syst.  Tech.  J.  50, 843-659  (1971). 

28ee,  also,  D.  Marcuse,  Theory  of  Dielectric  Optical  Waveguides, 
(Academic,  New  York,  1974). 

•’B.  Crosignani,  P.  Di  Porto,  and  C.  H.  Papas,  “Theory  of  time-de¬ 
pendent  propagation  in  multimode  lightguides,”  J.  Opt.  Soc.  Am. 
67. 1300-1307  (1977). 

4B.  Crosignani  and  P.  Di  Porto,  “Coherence  of  an  electromagnetic  field 
propagating  in  a  weakly  guiding  fiber,”  J.  Appl.  Phys.  44, 4616-4617 
(1973). 
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The  dependence  of  the  temporal  width  of  the  impulse  response  on  the  length  z  of  a  two-mode 
optical  fiber  is  examined.  This  quantity,  which  is  proportional  to  r  in  the  absence  of  mode  coupling 
end  to  z'ri  in  the  presence  of  weak  random  coupling  among  the  guided  modes,  possesses  a  different 
dependence  in  the  case  of  a  deterministic  resonant-coupling  model,  appropriate  for  describing  a 
rather  general  class  of  actual  situations.  The  relevant  role  played  by  the  coherence  time  of  the  signal 
is  demonstrated. 


INTRODUCTION 

The  influence  on  modal  dispersion  of  coupling  among  the 
guided  modes  of  a  multimode  optical  fiber  has  been  examined 
by  many  authors1  since  a  possible  reduction  of  the  pulse 
spreading  due  to  this  mechanism  was  first  predicted.2  The 
analysis  of  electromagnetic  propagation  in  a  fiber  in  the 
presence  of  mode  coupling  is  usually  accomplished  in  the 
frame  of  a  statistical  theory  which  deals  with  quantities  such 
as  the  powers  P„  of  the  individual  modes  averaged  over  a 
suitable  ensemble  of  similar  fibers.3  This  approach  allows 
one  to  investigate  the  electromagnetic  propagation  by  means 
of  a  set  of  differential  coupled  equations  describing  the  evo¬ 
lution  of  the  averaged  quantities  <P„ ),  which,  under  certain 
hypotheses,  can  be  cast  in  an  exactly  solvable  form.4  One  of 
the  main  results  of  this  theory  concerns  modal  dispersion. 
More  precisely,  it  turns  out  that  the  width  of  the  impulse  re¬ 
sponse  exhibits,  for  a  large  traveled  distance  z,  a  square- 
root- type  dependence  on  r ,  while  the  width  depends  linearly 
on  *  in  the  uncoupled  case.3  This  result  implies  that  for  long 
traveled  distances  mode  coupling  reduces  dispersion.  Re¬ 
duction  of  dispersion  due  to  mode  coupling  has  actually  been 
observed* 

The  statistical  approach,  which  has  some  intrinsic  limita¬ 
tions  due  to  its  very  nature,7  is  suitable  foe  dealing  with 
weak -coupling  situations  such  as  those  arising  from  the 
presence  of  many  small  unpredicub',  :.v. perfections  in  the 
fiber  structure.*  However,  there  are  libers  in  which  mode 
coupling  takes  place  in  such  a  way  that  they  cannot  be  con¬ 
sidered  as  typical  samples  of  the  ensemble  to  which  the  sta¬ 


tistical  theory  applies.  This  occurs,  in  particular,  in  con¬ 
nection  with  spatially  quasiperiodic  imperfections  giving  rise 
to  a  resonant  coupling  of  the  various  modes.  In  this  paper, 
we  wish  to  examine  electromagnetic  propagation  in  a  two¬ 
mode  optical  fiber  possessing  this  kind  of  coupling,  our  main 
aim  being  the  investigation  of  the  dependence  of  pulse 
spreading  on  length  traveled.  This  spreading  turns  out  to 
exhibit  a  behavior  that  is  different  from  the  square-root-type 
dependence  predicted  by  statistical  theory. 

The  two-mode  optical  fiber,  beyond  furnishing  the  simplest 
model  for  treating  mode  coupling,  is  interesting  per  se  due  to 
the  recent  developments  of  single-mode  optical  waveguides 
in  which  two  effective  modes  with  orthogonal  polarization  are 
present.910 


I.  PROPAGATION  IN  THE  PRESENCE  OF 
RESONANT  COUPLING  WITH  CONSTANT 
AMPLITUDE 

The  propagation  of  a  finite-bandwidth  field 
E(r,*,f)  ■  Ei(r)  J cl(z,u)e‘“‘~ie'M*dw 

+  Ei(r)  J" C2(r,«)e‘"<-'**<,'>,dw 

■  E,(r)o,  (x,t)  +  E2(r)o2(r,t)  (1) 
in  a  lossless  fiber  supporting  two  guided  coupled  modes  is 
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governed  by  the  system  of  equations3’7 

'  -j-cHz.u)  *  Ki2(z)eiM“uc2(z.w) 

dz 

■j-c2(z,w)  *  -lflt(*)e-iA<“)*ci(2,w), 


which  describes  the  evolution  of  the  mode  amplitudes  c,(z,w). 
We  have  indicated  with  r  and  z  the  transverse  and  longitu¬ 
dinal  coordinates,  with  A(w)  «  0i(w)  -  j32(w)  the  difference 
between  the  propagation  constants  of  the  two  unperturbed 
modes  Ei(r)  and  E2(r),  and  with  K\2(z)  the  coupling  coeffi¬ 
cient,  which  is  related  to  the  departure  of  the  fiber  structure 
from  the  ideal  configuration. 

If  we  assume  Ku(z)  to  be  of  the  form 

Ki2(z)  *  2Kcos(xr),  (3) 

with 

X  ■  A(w0),  (4) 

&>o  representing  a  frequency  within  the  spectral  width  5w  of 
the  signal,  the  set  of  Eqs.  (2)  becomes 

•—  Ci(z,w)  **  Xe‘lA<">+A<“»'l*C2(z,o)) 
dz 

+  Ke‘lA(“>-A<“»>lIC2(z.w),  (5) 

~  C2(z,w)  ■  -ff*e-*lA<“»-A<“®)l,Ci(z,4j) 
dz 

—  K*e-,lA(“’)+A(“o)llci(z,«). 

By  taking  advantage  of  the  rotating  wave  approximation 
which,  whenever  \K |  «  |  A(w)| ,  allows  us  to  neglect  the  rap¬ 
idly  oscillating  terms  of  the  kind  exp|±t[A(w)  f  A(oo0)]z|  that 
average  out  to  zero  over  the  distance  of  propagation,  the  set 
of  Eqs.  (5)  reduces  to 

d 

—  ci(z,«)  ■  Ke‘yM‘C2(z,u’), 

d  <6) 

—  c2(z,«)  -  -K*e-‘yM*ci(z,w), 
dz 

with  7(«)  “  A(w)  -  A(ojo).  The  set  of  Eqs.  (6)  can  be  easily 
solved  and  the  resulting  expressions  for  ci(z,w)  and  c2(z,w) 
inserted  into  Eq.  (1),  thus  giving11 

ai(z,t)  ■  exp(-iA(wo)z/2) 

X  J*  expj-iOriw)  +  /32(w)]z/2  +  iwt| 

X  U  exp(i|K|z  + »(w  -  wo)2z/8|K|02) 

+  B  exp(-i|K|z  -  i(w  -  «.0):z/8|K  |02]|dw  (7) 

and 

«a(*,t)  ■ » exp(-ip  +  »A(wo)z/2] 

X  J*  exp)-i(0i(w)  +  02iw)jz/2  +  iw£| 

X  \A  exp(«|K|z  +  Uu  -  wo)2z/8|K|02] 

-  B  exp|-»|K|z  -  »(w  -  wo)2z/8|K|02]|dw.  (8) 


Here, 

A  ■  (l/2)[ci(0,w)  —  te'*c2(0,w)]( 

B  ■  (l/2)[ci(0,u)>  +  ie‘*C2(0,«)l, 

where  exp(i$)  ■  K/\K\  and 

1/0  =  1 IV x-  MV 2, 

V}  *  (dfii/d(t>)umuo  indicating  the  inverse  of  the  group  velocity 
of  the  tth  unperturbed  mode.  Furthermore,  we  have  as¬ 
sumed 

£<i>«|K|0  (11) 

and 

A(io)  sat  A(w0)  +  (<*>-  w0 )/0,  (12) 

but  being  the  source  bandwidth. 

The  electromagnetic  powers  P\(z,t)  and  P2(z,t)  carried  by 
modes  1  and  2,  whose  behavior  determines  the  pulse  shape, 
are  proportional  to  (|ai(z,f)|2).vand  (|a2(z,t)|2>lv,  where  the 
symbol  (  },v  indicates  the  averaging  operation  over  the 
fluctuations  of  the  Bource  exciting  the  fiber.  In  order  to  de¬ 
termine  (|ai(z,t)|2),v  and  (|a2(z,£)|2>,v.  it  is  necessary  to 
evaluate  the  quantities  (ci(0,w)cl(0, <*/)>,„  <C2(0,w)cJ(0,w')}«v, 
and  (ci(0,a>)c2(0,u>')}«v  as  functions  of  the  source  bandwidth 
and  pulse  duration.  By  assuming  a  pulse  of  Gaussian  shape 
and  width  Tp  injected  into  the  fiber  by  a  source  possessing  a 
coherence  time  Tc  *  2ir/bu>,  one  has11 


(c,(0,u))c;(0,u)')>«v  -  (Ci(0,o)o)c1(0,a>')).v 
X  exp(— H27’i/4  -  Q'27?/4  +  QSl'Tl/2)  s  F(il,Q'),  (13) 


where 

7?  -  TUT j  +  T*)/(2T\  +  7?) 

(14) 

and 

7§  -  T*/(2T*  +  Tc), 

(15) 

with  fi  *  u 

-  w0  and  Q'  “  w'  -  w0. 

With  the  help  of  Eqs.  (7),  (8),  and  (13)  we  can  write,  after 
assuming  for  the  sake  of  simplicity  that  ci(0,w)  *  c2(0,w), 


<  |o,(z,t)|2).v  -  f  dfl 

X  H{z,Q,Q';K)F(W),  (16) 

where 

H(z,Q,Q';K)  -  <1/2)((1  +  sin*)e'<t,l-n'J>',<*> 

-Ml  -  sin*)e1<n',-°’>’’<*> 

+  i  cosde-2,'l*lt 

-  i  COS^e2‘lKl'  ei(Q,+  0'»)r»(*>j  (17) 

and 

r2(z)-z/8|tf|02  (18) 

the  same  expression  being  valid  for  ( |a2(z,£)|2>.v  but  for  the 
substitution  cos#  -*  -cos*. 

By  performing  the  double  integral  appearing  in  Eq.  (16), 
one  finally  obtains 


(9) 

1 

GO) 
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,lri  /T,irv  _  4t<Ci(0,wo)c!(0,wo)).v  .  „  „/fr,  .  ..  4»cos^e,(0,uo)Ci(0,wo)).v 

< K*.'>l  >.v  nm  +  16r  w|)l/2 «Pl"2‘ 2/^  +  16t'<*,/7^  *  fiTifTTeHi,??^^ 

X  exp(-2t'2(7in  +  167tT*(2)]/({7lTg  -  16r<(z)]2  +  647V(z)|) 

X  sin(2t'*|47VU)  +  64t«U)J/|[TST?  -  16r<(z)]2  +  647V(*)I  +  2|/C|z  +  */2),  (19) 


4*  cos(»(cl(0,h)o)ci(0,mo))< 


where 

t'  -  t  -  z/V,  (20) 

with 

1/V  -  (1/Vi  +  l/V2)/2,  (21) 

tl~Til-Ti’-Tffi/(2T2p+Tl:),  (22) 

r*»7t+7l=7%  (23) 

and 

t  -  arg[7lTg  -  16r*(z)  -  8»T?r2(z)].  (24) 

The  significant  result  contained  in  Eq.  (19)  is  that  the 
powers  P\(z,t )  and  P2(z,t )  exhibit  the  same  space-time  evo¬ 
lution,  traveling  with  a  bulk  velocity  V  and  undergoing  a 
broadening  which,  for  z  sufficiently  large  (z  >  z,  z  being  de¬ 
fined  through  the  relation  r2(z)  *  TpTJ 4),  is  proportional 
to  the  traveled  length  itself. 

Actually,  it  is  the  first  term  on  the  right-hand  side  of  Eq. 
(19)  that  plays  the  leading  role  for  what  concerns  broadening, 
since  T3  <  Tp  (see  Eqs.  (22)  and  (23)].  In  practice,  Te  «  Tp. 
so  that  T 3  at  TJyfi  «  Tp.  Accordingly,  the  temporal  width 
Tp(z)  of  the  pulse  can  be  evaluated  as 

Tp(z)  »  z/(2|ff|i52Tc)1/2  for  z  >  z,  (25) 

so  that  the  pulse  temporal  width  behaves  linearly  with  z,  as 
in  the  absence  of  mode  coupling. 

The  width  Tp(z)  can  be  easily  shown  to  be  much  smaller 
than  T  “  z/0  relative  to  the  noncoupling  case,  that  is,  to  the 
difference  between  the  group  delays  of  the  two  unperturbed 
modes.  As  a  matter  of  fact,  one  has 

Tp(t)  -  (77v'2)(B«/2t|K|i»,  (26) 

so  that,  according  to  Eq.  (11),  Tp(z)«T.  This  fact  confirms 
the  positive  role  played  by  coupling  in  reducing  dispersion. 
Thus,  the  influence  of  mode  coupling  on  pulse  broadening  is 
only  of  a  quantitative  nature  in  our  deterministic  model,  while, 
in  the  framework  of  the  statistical  theory,  mode  coupling  re¬ 
sults  in  a  z  dependence  of  the  pulse  width  of  the  square-root 
type.  This  latter  kind  of  dependence  can  be  accounted  for 
in  an  intuitive  way  by  taking  advantage  of  the  Central  Limit 
Theorem  in  connection  with  the  random  “jumps"  of  each 
photon  from  one  mode  to  the  othtr.'2  «•!  n<>  random  process 
occurs  in  our  model. 

We  wish  to  conclude  this  sectioi.  b>  noting  the  relevance  of 
the  coherence  time  Te  of  the  source  (provided  by  Eq.  (11)  is 
verified]  in  determining  the  pulse  spreading  in  our  case,  while 
T,  (or,  equivalently,  the  source  bandwidth)  does  not  influence 
the  pulse  shape  according  to  the  rm.-I-  ns  of  the  rtatistical 
approach.3  This  is  directly  related  to  the  basic  hypothesis 
underlying  the  statistical  theory,  according  to  which  the  field 
amplitudes  do  not  vary  appreciably  over  the  small  distance 
after  which  the  mode-coupling  random  imperfections  decor¬ 
relate,  leading  to  a  system  of  coupled  equations  involving  only 


the  mode  powers  <  P„ ) . 3  On  the  contrary,  coherence  effects 
are  preserved  in  our  situation  due  to  the  periodicity  of  the 
coupling  coefficient,  a  fact  that  does  not  allow  consideration 
of  our  fiber  as  a  typical  sample  of  the  statistical  ensemble. 


II.  INFLUENCE  OF  VARIATIONS  OF  COUPLING 
AMPLITUDE  AND  PERIODICITY 

A  more  general  coupling  model  can  be  obtained  by  letting 
the  amplitude  K  and  the  periodicity  x  of  K12(z)  to  assume 
different  constant  values  K,  and  Xi  over  different  fiber  regions 
(possessing,  in  general,  different  lengths  L,).  In  order  to 
preserve  the  resonance  condition,  we  assume  the  existence  of 
frequencies  Ui  lying  inside  the  bandwidth  ioj  such  that 

Xi  -  0i(w;)  -  02(«,).  (27) 

From  an  analytical  point  of  view,  the  problem  of  pulse 
propagation  in  this  structure  can  be  easily  solved  by  observing 
that,  due  to  the  form  of  Eq.  (16),  the  quantity 

H(Ll,Q,n';Ka)F(Q,Q')  (28) 

furnishes  the  “boundary  condition”  for  deriving  the  solution 
in  the  second  region,  which  reads 

<|a,(z,t)|2>.v« 

X  J/(z2,fl,fl,;K2)/f(L),0,Q';if1)F(fl,fl,)1  (29) 
for  0  £  z2  £  Li  and  z  *  L\  +  z2.  More  in  general,  one  has 

<|aj(z,f)|2>.v  - 

X  H(zn,Sl,W;Kn)H(Ln-i,Q,(l';Kn-\) 

X  H{L„-2,W;Kn-2)  X  •  ■ .  H(Ll,fl,n';K1)F(Q,fl').  (30) 

for  0  S  z„  £  L„  and  z  *  Li  +  £,2  +  •  •  •  +  L„-i  +  *n- 

The  right-hand  side  of  Eq.  (30)  is  the  sum  of  a  number  of 
terms  of  the  kind 

J*  dftj*  dfl'exptiafl2  +  tbfl'2) 

X  exp[i(fl  -  »')«  -  z/V)]F(H,Q'),  (31) 

where  a  and  b  are  linear  combinations  of  the  quantities  r2(Li), 
r2(L2), . . .,  t2(z„)  with  coefficients  +1  and  -1  in  all  possible 
combinations.  By  recalling  Eq.  (13),  the  double  integral  ap¬ 
pearing  in  Eq.  (31)  can  be  explicitly  evaluated.  The  most 
relevant  contributions  to  ( |ai(z,f)|2>,v  come  from  the  terms 
for  which 

fl  -  -b  -  ±(r2(L,)  +  rHLt)  +  •  -  *  +  r2(z„»,  (32) 

which  evolve  in  the  form  of  a  pulse  moving  with  a  velocity  V 
and  possessing  a  temporal  width  Tp(z)  »  TP(L\  +  L%  +  •  *  • 
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+  Lh-i  +Zn)  given  by 
TPU) 

Li/|K,1  +  Lj/IKal  +  •  ■  •  +  L„-i/|Kn-,|  +  zJ\Kn 
V20*Te 


Equation  (33)  generalizes  Eq.  (25)  but  assures  a  linear  de¬ 
pendence  of  TPU)  on  the  traveled  length  only  inside  each 
single  fiber  region  with  a  given  value  of  |  K  | . 


IU.  CONCLUSIONS 

We  have  investigated  the  broadening  of  a  pulse  propagating 
in  a  two-mode  fiber,  made  up  of  a  number  of  regions  in  each 
of  which  a  spatially  periodic  resonant  coupling  with  constant 
amplitude  is  present  This  kind  of  optical  waveguide  does  not 
belong  to  the  statistical  ensemble  considered  in  the  random¬ 
coupling  theory,  so  that  electromagnetic  propagation  cannot 
be  treated  by  means  of  coupled-power  equations,  and  the 
mode-amplitude  evolution  has  to  be  explicitly  investigated. 
As  a  consequence,  the  temporal  broadening  of  the  pulse  turns 
out  to  be  related  to  the  source  coherence  time,  while  it  does 
not  exhibit  the  dependence  on  the  square  root  of  the  fiber 
length  predicted  by  the  statistical  approach. 
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SUMMARY 


V. 
\ 

\ 


The  propagation  in  a  multimode  optical  fiber  of  a  finite-bandwidth  optical  carrier 
aodulated  by  a  nonstationary  signal  is  investigated.  The  fluctuations  of  the  field  due  to 
random  node-couoling  are  considered  and  the  set  of  counled  equations  describing  their  evo_ 
lution  is  derived.  In  particular,  this  allows  us  to  investigate  the  propagation  of  a  fre_ 
quency-modulated  signal  and  to  obtain  a  general  theorem  concerning  the  asymptotic  behavior 
of  mode-power  fluctuations. 

1 .  INTRODUCTION 


The  propagation  of  a  stationary  monochromatic  carrier  in  a  multimode  optical  fiber 
in  the  presence  of  mode-coupling  has  been  extensively  studied  by  Marcuse  (1),  in  the  fra_ 
ms  of  a  statistical  approach  based  on  the  introduction  of  an  ensemble  of  similar  fibers 
slightly  differing  from  a  common  ideal  structure  for  the  presence  of  random  imperfections. 
The  procedure  allows  one  to  write  two  closed  systems  of  differential  coupled  equations  for 
the  single-mode  powers  averaged  over  the  ensemble  of  fibers  <  P  y  and  for  the  correlati_ 
ons  <  P  P  >  (0,(2).  The  extension  to  the  nonstationary  nonamnochromatlc  ease  has  been 

accomplished,  for  what  concerns  the  average  powers  <P  >  ,  on  an  intuitive  basis  by  Mar_ 
cues  (1)  and  by  means  of  more  rigorous  treatments  by  Personlck  (3)  and  Steinberg  (4). 

We  generalise  the  statistical  approach  to  furnish  the  evolution  of  the  quantities 

<e  (s,«0)c"(s,  >  and  <  c  (s,  iO)c*(z,  bi  ’)  c  (s,  u)")e»(s,  a)  •")  >  ,  where  the  c  *s 

■  r  m  n  r  K  ra 

are  the  slowly  varying  mode-amplitudes  and  the  bar  denotes  the  averaging  operation  over 

the  souroe  fluctuations.  The  knowledge  of  these  quantities  allows  us  in  turn  to  describe 

the  behaviour  of  the  field  and  of  its  fluctuations  in  the  general  situation  in  which  a  fi_ 

nite- bandwidth  carrier  is  amplitude  or  frequency-modulated. 

An  other  significant  result  concerns  the  asymptotic  behavior  of  the  single-mode  po_ 
w*r  fluctuations,  which  tend  to  vanish  as  a  consequence  of  the  finite  carrier-bandwidth, 
so  that,  in  the  lossless  case,  power  eouipartition  among  modes  takes  place  not  only  on  the 
average,  but  over  each  fiber  of  the  ensemble  (5). 

t 

I  2.  DESCRIPTION  OF  ELECTROMAGNET  I C  PROPAGATION  IW  A  FIBER 

f 


I  The  transverse  part  of  the  electric  and  magnetic  fields  excited  in  a  fiber  by  means 

of  a  finite-bandwidth  source  can  be  approximately  written  in  terms  of  the  forward  traveling 
|  guided  modes  E^  ,  as  (1) 


dence  of  the  slowly-varying  mode  amplitudes  cr  la  associated  with  the  departure  of  the  ft, 
her  structure  from  the  ideal  one  to  which  the  configurations  E  ,  R  pertain.  The  eleotro_' 
magnetic  cower  carried  through  an  arbitrary  portion  of  area  O’- ™ ofTny  given  fiber  section 
a  »  const,  can  be  expressed  in  terms  of  the  complex  Poynting  vector  S  (6) 

s  -  i  Ex  H*  ,  (4) 
where  the  bar  Indicates  the  averaging  operation  over  the  tuu.ue  fluctuations,  as 


,r  ■  8e  »  .  S  dr 

J 


(5) 


Here  "Re"  means  "real  part  of"  and  e  is  the  unit  vector  in  the  positive  s-dtreo_ 
tion.  By  means  of  Eos.(l),(2)  and  (4),  Eq.75)  yields  “ 


where 

A  F  - 
mr 


prU.t)  «  *X  X 

ra  r 

(  (  i(d  -  (J'H-iffi  (« )- A  (W)ls _ 
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The  total  power  P  carried  by  the  field  is  obtained  by  integrating  e  .  S  in  Ea.(5) 
over  the  whole  plane  a  »  const.  .  This  operation  furnishes,  with  the  help- Sf~ the  ortho_ 
gonality  relation  between  the  modes  (1)  ” 

> oo  f  too 


jdx  yy  vVS'^o^^'^o'  -  2  pSmn 

-  o»  -oo  r 

whers  P  is  a  positive  normalization  constant  and  A  la  the  usual  Kronecker  symbol, 

nra 
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where 
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ID 


-  i[&  (fa>)  -(S  (W)l*  _ , _ 
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(10) 


can  be  interpreted  as  the  cower  carried  by  the  nth  mode.  It  is  worthwhile  to  note  that 
the  interference  terms  between  the  various  modes  disappear  in  the  expression  of  the  to 
tal  power  P^,  while  they  are  nresent  in  that  of  P  •“  . 

3-  STATISTICAL  DESCRIPTION  OP  PROPAGATION 


The  random  nature  of  the  imperfections  unavoidably  present  in  a  real  fiber  suggests 
some  statistical  procedure  aa  the  most  natural  way  for  describing  the  propagation  in  opti_ 
cal  fibers.  This  is  accomplished  by  introducing  an  ensemble  of  similar  fibers,  each  of 
Which  differs  from  a  common  ideal  model  for  the  presence  of  random  imperfections,  and  to 
evaluate  the  average  of  the  significant  ouantities,  e.g.  P*”,  P  and  P_,  over  this  ensem 

^ ^  at  *“ 

ble,  whtch  In  turn  is  equivalent  to  consider  <  c  (z,  h>  £•)')■>  as  the  basic  quantity. 

!fl  r 

In  order  to  test  the  relevance  of  this  approach  for  what  concerns  practical  aitua_ 
tions,  in  which  a  single  fiber  Is  usually  employed,  one  has  to  investigate  the  statisti_ 
cal  fluctuations  around  the  average  values.  This  leads  to  consider  higher  order  averages 

of  the  kind  <  c  (z,  )c*(  z,  tO')  c  ( i ,  cJ  "  >  c*  (  z ,  «•"  )>,  <  c  ( z  c"(  t  ,Q' )  c  (z,fa)")c*(z,tJ"'  )>. 

in  r  n  r  ^  id  r  r  n 

One  starts  from  the  system  of  equations  describing  the  evolution  of  the  c  's  in  the 

single  fiber,  which  reads  (1)  m 


de  («,u>) 
w 

da 


22. 


with 


A^z.o))  ck(*,W) 


(11) 


— >  ' 


Amk(*,W  ) 


K«k(B) 


#i[  (&.(*>)  -  (ik(tf )]  * 


(12) 


the  V’8  being  coupling  coefficients,  vanishing  in  the  case  of  an  ideal  fiber. 

By  using  Bn. (11)  one  can  easily  derive  the  following  equations 
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+  Amk*B,W’^  cn(z»W)ck(*.tt)')  }  .  (13) 
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ediieh  have  to  be  averaged  over  the  fiber  ensemble.  The  resulting  enuations  can  be  out  in 
the  form  of  closed  systems  under  the  homogeneity  assumntlon 


<  «■,<*>  >  •  0  .  (15) 

<  V*r.<* <17 1  '  •  ('«) 

provided  that  the  counting  is  small  enough  that  the  c  ’s  do  not  significantly  vary  over 
the  correlation  length  of  the  statistical  variables  K®  .  In  this  way  one  obtains  (see  also 
Ref. (7)  ),  as  a  first  significant  result,  an 

3z  <  ca(*-fci1«;(*.  «’)>--  ^  <  eB(z,<0  )c»(z,u>’)  >  ,  m/r  ,  (17) 
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The  real  part  of  g^  can  be  shown  to  be  a  positive  quantity  (7),  so  that 
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vanish  for  i  large  enough  (  or  are  identically  aero  if  they  are  such  at  a-O). 

The  lowest-order  eouations,  concerning  single-mode  powers,  read  (see  also  Refs. (3) 
and  (4)  ) 


fj  <«„<«.«  )c;<«.tf*>> 


-«Z 

v« 


h^) -C  cB(8,«)c*(a,aj')> 


V  t  F  tJ)-  ^VC  60*>-  >3  a  - 

♦  /  h.e  I*  rm  ’  k  rit  ^  ek(a,w)cj(«,4Ji 

kj& 


>  ,  (23) 


c°° 

J<K«1c7)K1 


ifflW  - 


which  constitutes  a  closed  set  of  eauations. 

The  higher-order  terms,  concerning  mode-power  fluctuations  and  mode-mode  correla_ 
tions,  cannot  he  considered  separately  since  they  form  together  a  closed  set  of  equations, 
which  reads 
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4.  SOME  APPLICATIONS 


The  set  of  equations  governing  the  behavior  of  the  average  mode-power  <  P  (z,t) > 
can  be  easily  obtained  by  taking  advantage  of  Eos. (10)  and  (23),  thus  getting  (3)?(4),(7) 


-~<PB(z,t)>  +  ^  7§t  <V‘,t);>a^  hak  fcPk(«,t>>  -  <Pa(z,t)>]  , 
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where 


(27) 


(28) 


is  the  group  velocity  of  the  mth  mode.  On  the  other  hand,  Eqs.(6),(7)  and  (10)  allow  one 
to  write 
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where  the  <  Pm(z,t)  >  possess  the  same  soace-tlae  dependence  as  the  <P  (z,t)  >  .  Sin_ 
ee  a  frequency  modulation  of  the  signal  does  not  affect  the  ■<  P  >  's,  th?  dependence  on” 
this  modulation  of  <  P*"( z , t )  >•  is  contained  in  the  <  A®J>  ’s, "which  vanish  when  consi_ 
dering  the  power  transmitted  through  the  whole  section  or  the  fiber.  Prom  Eqa. (7)  and  (17) 
it  follows  that  (7) 
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where 

-i<0  t 

P(t)  -  E( z*0,t)  e  0  .  (31) 

Equation  (30)  shows  that  mode-coupling  affects  the  propagation  of  a  frequency-mo_ 
dulated  signal  in  an  ensemble  of  fibers  only  for  the  presence  of  a  spatial  damping  factor 
characteristic  of  each  couple  of  modes.  More  precisely,  if  one  excites  at  the  fiber  in_ 
put  only  two  modes  a  and  r,  all  interference  terms  but  <  >  vanish  throughout  the  fT_ 

her  (see  Eo.(17)),  while  the  surviving  term  propagates  as  vh  the  absence  of  coupling  but” 
for  the  spatial  attenuation. 

The  set  of  Eos. f 27 ^  and  the  set  of  Eos. (25)  and  (26)  allow  one  in  principle  to  de 
t ermine  the  quantities  <  P  (z,t)>  ,  £  P*(z,t)  >  and  <P  (s,t)P  (s,t)  >  ,  and  thus  the 
fluctuations  whose  magnitude  furnishes  a  criterion  for  the "applicability  of  the  statisti_ 
cal  theory  to  a  single  fiber. While  the  general  solution  of  Eps. (25)  and  (26)  is  a  formidi 
ble  task,  an  Important  particular  result  can  be  obtained.  If  one  Introduces  the  energy 
I  (s)  pertaining  to  the  mtr.  mode 

IB  ^  ^ 

y*)  *  j  PB(a,t)dt  -  2TTP  J  cB(z,«)c;(t,fc))  d  (A 

"o* 

It  can  be  shown  (7),  by  means  of  Eqs.(25),  that,  under  the  condition 


(32) 
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and  provided  that  z  exceeds  the  soul partition  distance  beyond  which  <■  1^  > 
the  following  relation  holds 


<*l> 

<  I->2 


<<  1 


<ln>' 


(34) 


which  ladles  that  the  statistical  uncertainty  of  the  ou entity  I  (a)  over  the  fiber  en_ 
seable  is  negligible,  so  that  aode-energy  ecu! partition  takes  place  over  each  fiber.  This 
last  consideration  holds  true  in  the  absence  of  losses  or,  more  in  general,  If  losses  are 
not  effective  over  the  characteristic  length  after  which  Eo. (34)  Is  verified. 
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